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THEORY  OF  IONIZING  ATOMIC  COLLISIONS 
AT  THERMAL  ENERGIES 

By 

LAWRENCE  RALPH  BELYEA 
June  1980 

Chairman:  D.  A.  Micha 
Major  Department:  Chemistry 

Ionization  in  low-energy  atomic  collisions  involving 
electronically  excited  species  is  investigated  using 
coupled  channels  and  effective  Hamiltonian  theories  for 
the  heavy  particle  motion.  Within  the  coupled  channels 
formulation  detailed  electron  energy  distributions  are 
obtained  for  the  Penning  ionization  of  Ar  by  He(1s2s,3S); 
the  structure  is  interpreted  as  arising  from  final  channel 
heavy  particle  resonances,  and  this  interpretation 
provides  new  insight  into  the  role  of  the  nuclear  dynamics 
on  electron  emission.  The  Riccati  transform  of  the  matrix 
SchrOdinger  equation  arising  in  the  coupled  channels 
theory  is  solved  to  obtain  cross  sections  for  Penning 
ionization  over  a range  of  energies,  and  results  were 


IX 


found  to  compare  well  with  experiment.  The  validity  of  a 
simplifying  decoupling  approximation  is  demonstrated  by 
numerical  studies;  the  agreement  between  "exact"  and 
decoupled  results  is  excellent  over  a range  of  energies, 
and  is  independent  of  heavy  particle  angular  momentum 
partial  wave. 

An  effective  Hamiltonian  theory  for  Penning  and 
associative  ionization  is  formulated  using  projection 
operator  techniques;  an  effective  potential  is  obtained 
for  the  nuclear  motion  of  a discrete  electronic  state 
which  is  unstable  with  respect  to  decay  by  electron 
emission.  Analysis  of  a complex,  non-local  and 
energy-dependent  potential  energy  kernel  demonstrates  for 
the  first  time  that  while  a local  approximation  to  the 
real  part  may  be  valid,  such  approximation  may  not  be 
valid  for  the  imaginary  part. 

A method  is  presented  for  solving  inhomogeneous, 
second-order  ordinary  differential  equations  with  a 
complex,  non-local  kernel,  and  simplifications  for  the 
case  of  a purely  imaginary  kernel  are  discussed. 
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INTRODUCTION 


The  interaction  of  electronically  excited  atoms  and 
molecules  with  ground  state  species  is  of  practical 
interest  in  lasers,  plasmas  and  atmospheric  chemistry  and 
of  theoretical  interest  for  the  understanding  of  energy 
transfer  mechanisms.  The  subject  of  this  dissertation  is 
the  quantum  mechanical  study  of  the  dynamics  and 
energetics  of  ionizing  collisions  in  atoms  where  the 
nature  of  the  energy  transfer  is  in  the  form  of  electronic 
excitation;  in  particular,  attention  is  focused  on  the 
ionization  of  a neutral  target  atom  in  its  ground 
electronic  state  by  an  electronically  excited,  metastable 
atom. 

In  addition  to  the  elastic  scattering,  various 
reactive  channels  or  pathways  are  available  depending  upon 
the  collision  energy,  the  excitation  energy  of  the 
metastable  atcm,  and  those  excitations  and  ionizations  of 
the  target  which  are  allowed  by  energy  conservation 
considerations.  If  the  metastable  atom  excitation  energy 
exceeds  the  ionization  potential  of  the  target,  ionization 
can  occur;  if  this  excitation  is  less  than  the  target 
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ionization  energy,  excitation  and  ionization  of  the  target 
are  both  possible  modes  of  energy  transfer.  Final 
channels  (or  post-collision  rearrangements)  will  in 
general  then  consist  of:  1)  elastically  scattered 
reactants,  2)  de-excited  metastafcle  and  excited  target,  3) 
de-excited  metastable,  ionized  target  plus  electron,  and 
4)  molecular  ion  plus  electron. 

It  is  the  object  of  this  work  to  put  forth  a theory 
of  ionizing  collisions  of  atoms  with  metastables  in  the 
case  where  the  electronic  excitation  energy  of  the  excited 
atom  exceeds  the  ionization  energy  of  the  target;  the 
processes  are  referred  to  as  associative  or  Penning 
ionization  (AI  or  PI)  depending  upon  whether  the  heavy 
particles  are  bound  or  unbound,  respectively,  in  the 
ionized  channel.  Such  a theory  should  be  sufficiently 
detailed  so  as  to  be  able  to  be  predictive  of  some  of  the 
wide  range  of  experimental  observations  which  are  possible 
for  such  a collision  process:  roughly  in  order  of 
increasing  detail,  one  may  measure  total  ionization  cross 
section  and  its  energy  dependence,  branching  (AI/PI)  ratio 
and  its  energy  dependence,  angular  distribution  of 
elastically  scattered  particles,  angular  distribution  of 
de-excited  metastaole  and  ionized  target,  vibrational  and 
rotational  distributions  for  AI,  and  angle-  and  energy- 
dependent  distributions  of  the  ejected  electron.  Analysis 
of  results  such  as  these  can  provide  useful  insights  into 
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features  of  both  the  collision  dynamics  and  the  ionization 
raecha  nism. 

Following  a mathematical  review,  a brief  survey  of 
approaches  to  the  collisional  ionization  problem  will  be 
made  in  Chapter  I.  To  complete  the  review  and  anticipate 
some  of  the  numerical  results  to  be  presented  later,  a 
recent  coupled-channels  approach  which  makes  use  of  a 
discrete  approximation  to  the  electronic  continuum  will  be 
discussed. 

In  Chapter  II  a discussion  of  the  electronic  states 
appropriate  to  the  problem  of  ionizing  heavy-particle 
collisions  will  be  presented.  The  construction  of  proper 
molecular  resonance  and  continuum  electronic  states  is 
discussed  in  some  detail.  Then  the  collision  problem  for 
nuclei  and  electrons  will  be  analyzed  in  terras  of  a 
projection  operator  formalism;  given  a basis  of 
Born-Oppenhei mer  electronic  states,  the  space  spanned  by 
these  states  is  partitioned  in  a physically  appealing  way 
into  a "Q-space"  or  "bound"  portion  and  a "P-space"  or 
"continuum  (scattering)"  subspace.  Heavy  particle  motion 
in  the  entrance  channel  (whose  state  has  a zero  projection 
in  P)  will  then  be  described  in  terras  of  a complex, 
non-local,  energy-dependent  effective  Hamiltonian  which 
accounts  for  both  elastic  scattering  and  ionization.  The 
various  transition  amplitudes  relevant  to  AI,  PI  and 
elastic  scattering  will  be  derived. 


along  with  the 
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corresponding  scattering  cross  sections.  The 
computational  scheme  developed  for  the  solution  of  the 
scattering  eguations  in  a coupled  channels  theory  will  be 
discussed  in  Chapter  III.  Application  of  the  theory  for 
Penning  ionization  will  be  made  for  a particular  choice  of 
collision  partners,  He*(1s2s,3S)  and  Ar,  and  various  cross 
sections  and  energy  distributions  will  be  compared  to 
recent  experimental  results.  Analysis  of  electron  energy 
distributions  is  seen  to  lead  to  a more  refined  picture  of 
the  Penning  "mechanism"  in  terms  of  a certain  type  of 
heavy-particle  resonance. 

In  Chapter  IV  we  will  discuss  the  solution  of  an 
inhomogeneous  differential  equation  describing  heavy 
particle  motion  for  Penning  and  associative  iorization  in 
atomic  systems.  Various  approximations  to  the  real  and 
imaginary  parts  of  the  complex,  non-local  kernel  are 
examined  in  detail.  We  will  present  a method  for  the 
solution  of  inhomogeneous  differential  equations  for  such 
a kernel;  the  treatment  is  general,  and  we  comment  on  its 
implementation  for  a particularly  useful  approximation  to 
the  kernel.  The  chapter  concludes  with  a summary  of  the 
present  work  and  a discussion  of  problems  where  further 
study  is  warranted. 


CHAPTER  I 
PRELIMINARIES 

1.1  Separation  of  Electronic  and  Nuclear  Motion 
Our  goal  lies  in  the  solution  of  the  Schrttdinger 
equat ion 

(E  - H)  V = 0 (1.1) 

where  H is  the  Hamiltonian  for  a system  of  nuclei  and 
electrons.  Inasmuch  as  this  work  is  devoted  to  a study  of 
ionizing  collisions  between  atoms,  we  restrict  our 
treatment  in  what  follows  to  diatomic  systems.  We  write, 
for  an  N-electron  system, 

H = H + T 

el  R (1.2) 

with  T , the  nuclear  kinetic  energy  operator,  given  by 
R 

T = ~(*2/2M)  V 2 

R R (1.3) 

M = HM/(M+M) 

A B A B (1.4) 

where  a (M  ) is  the  mass  of  atom  A (B) . The  electronic 
A B 

Hamiltonian  H is 
el 

H = Z(i=1,N)  {(-J<2/2m  ) V 2 

el  i i 

-Z  e2/ | r -r  | - Z e2/| r -r  | 

A i A B i B (1.5) 
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+ Z(j<i=1,N)  e2/ J r -r  |}+Z  Z e2/|R  - R |. 

i j A B A B 

The  terms  are  the  usual  electron  kinetic  energy  and 
electron-nuclear,  electron-electron  and  nuclear-nuclear 
interactions.  We  have  neglected  both  the  mass 

polarization  terms  which  arise  from  the  transformation  to 
center  of  mass  coordinates  (Froman,  1962)  and  spin  orbit 
interactions.  The  splitting  of  H (1.2)  is  done  to  exploit 
the  mass  disparity  between  electrons  and  nuclei,  in  order 
that  we  may  solve  separately  for  electronic  and  nuclear 
motion. 

We  now  introduce  a complete  and  orthonor mal  set  of 

discrete  electronic  states  [<J>  (r;P)}  where  r denotes  the 

i 

collection  of  electronic  coordinates,  R is  the 
internuclear  separation,  and  i labels  a particular  state. 
These  states  are  adiabatic  states  in  the  sense  that  their 
dependence  upon  the  internuclear  coordinate  is  parametric 
rather  than  dynamic.  Apart  from  this  restriction,  these 
functions  are  otherwise  unspecified  at  this  point.  We 
merely  mention  here  that  this  basis  must  be  augmented  by 
continuum  electronic  states  when  we  wish  to  study  the 
collisicnal  icnization  problem  in  the  next  chapter. 

Expanding  If  in  this  basis  (our  notation  is  that  of 

f(r,R)  - £(i=1,n)  <)>  (r ; R)  X (R) 

i i 


Tully  ( 1976)) 
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= Z(i=1,n)  <t>  X 

i i (1.6) 


substituting  into  (1.1),  premultiplying  by  4>  and 

j 

integrating  over  electronic  coordinates  yields,  for  the 

expansion  coefficients  (i.e.,  nuclear  wavef unctions) , 

(E  - T - V - T"  ) X = 

H jj  jj  j 

Z(i*j=1,n)  (V  +T"  +T ' )X 

ji  ji  ji  i (1.7) 

where 

T"  = - (X2/2M)  <$  |V  2<J>  > 

ji  j R i (1.8) 

T*  = - (X2/2M)  <<|>  | V <(>  >°l 

ji  j R i R (1.9) 

V = <$  |H  | <J>  > 

ji  j el  i (1.10) 

and  where  <. ..>  denotes  an  integration  over  electronic 
coordinates;  we  have  used  the  fact  that  for  ()>  real  and 


ortho normal,  the  diagonal  terms  T'  are  identically  zero 

jj 

(Dalgarno  and  McCarroll,  1956;  1957).  The  diagonal  terms 

V are  the  potential  energy  surfaces  for  nuclear  motion; 

ji 

the  terms  T"  are  non-adiabatic  corrections  to  these 

j j 


surfaces. 

Non-adiatat ic 

in  teractions 

are  given 

by 

the 

terms  T" 

and  T* 

All  T"  and 

T*  terms 

may 

be 

ji  ji 


collectively  thought  of  as  describing  the  response  of  the 
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electrons  or  of  the  electronic  state  to  the  nuclear 

motion.  The  result  (1.7)  is  exact,  given  the  form  (1.2) 

of  H and  the  adiabatic  hypothesis. 

The  mass  disparity  between  electrons  and  nuclei  has 

as  a consequence  a qualitative  difference  between 

electronic  and  nuclear  motions:  electronic  motion  is 

comparatively  rapid  and  highly  quantized,  while  the  heavy 

particle  motion  is  considerably  slower  and  more  nearly 

classical.  This  observation  is  the  essence  of  the 

Born-Oppenheimer  approximation,  which  mathematically 

amounts  to  neglect  of  all  T"  and  T'  terms  in  the  last 

result.  This  gives 

(E-V  -T  ) X = Z(i*j  = 1,n)  V X . 

jj  R j ji  i (1.11) 

Justification  for  this  step,  based  upon  order-of-ma gni tude 

estimates  of  the  quantities  involved,  may  be  found  in 

O'Malley  (1971).  Without  further  specification  of  the 

electronic  basis  we  see  that,  even  in  the  B-0 

approximation,  transitions  among  states  are  possible  by 

the  couplings  through  H . With  this  we  have  completed 

el 

our  review  of  the  separation  of  nuclear  and  electronic 

motions.  we  now  turn  our  attention  to  a discussion  of  the 

possible  choices  for  the  states  (4>  } . 

i 
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1.2  Representations  of  electronic  Wavef unctions 
We  have  as  yet  left  unspecified  the  set  of  electronic 
states  {$}  . There  are  two  classes  of  representations  which 
have  been  utilized  in  descriptions  of  atomic  collision 
processes;  these  are  the  adiabatic  representation  on  the 
one  hand,  and  a variety  of  diabatic  representations  on  the 
other. 


The  standard  adiabatic  representation  (Smith,  1969) 

of  quantum  chemistry  is  obtained  by  choosing  for  the  set 

{(}  those  states  which  are  eigenfunctions  of  H , i.e., 

el 

H 4>  = E <p  , 

el  i i i (1.12) 

where 

E (R)  = < 4 |rt  | <t>  > = V S 

i i el  i ij  ij  (1.13) 

The  result  (1.7)  now  reads 

(E  - T - E - 1"  )X  = 

R j j j 

£(i# j=1,n)  (T"  + T'  ) X . 

ji  ji  i (1.14) 

Couplings  through  H are  absent;  the  off-diagonal  T"  and 

el 

T'  terms  are  the  couplings  which  account  for  transitions 
to  other  energy  surfaces.  During  a collision  process  the 
particles  can  be  thought  to  move  along  the  resulting 
adiabatic  potential  energy  surfaces,  and  transitions  take 
place  between  states  due  to  the  nuclear  kinetic  energy 
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operator.  In  regions  where  adiabatic  surfaces  are  in 
close  proximity,  the  otherwise-small  couplings  can  become 
quite  large  and  it  is  in  these  regions  that  there  is  a 
high  probability  of  transition  from  one  state  to  another. 
If  we  invoke  the  B-0  approximation,  we  obtain,  in  analogy 
with  (1.11), 

(E-T  - E ) X = 0. 

R j j (1-15) 

At  this  point  the  states  X have  become  completely 

j 

uncoupled,  so  what  we  have  is  a representation 
corresponding  to  stationary  states  of  molecules.  The 
states  of  the  standard  adiabatic  representation  exhibit 
another  important  property,  namely  that  the  non-crossing 
rule  (von  Neumann  and  Wigner,  1929)  applies  to  states  of 
the  same  symmetry.  This  results  in  "avoided  crossings"  of 
potential  energy  surfaces,  together  with  large  electronic 
wavefunction  gradients  in  the  vicinity  of  these  points. 
Although  this  description  has  enjoyed  widespread  use  by 
quantum  chemists,  it  is  clear  that  it  is  wholly  unsuitable 
for  the  scattering  theorist  when  considering  a reactive 
collision  problem  where  the  system  can  no  longer  be 
thought  to  stay  on  one  PE  surface.  It  is  conceivable  to 
calculate  adiabatic  energy  surfaces  and  all  the  B-0 
coupling  terms  using  electronic  eigenstates  and  then 
transform  to  a different  representation  more  suited  to  a 
particular  collision  problem.  However,  electronic 


wavef unctions  are  usually  determined  for  a few  select 
clamped-nuclei  geometries;  since  these  wavef unctions  do 
not  have  explicit  dependence  upon  the  internuclear 
coordinates,  and  as  these  coordinates  are  actually 
continuous  variables,  it  is  easy  to  see  that  this  is  not  a 
feasible  program  to  carry  out.  An  added  consideration  is 
the  size  of  the  expansion  basis  {0}  and  hence  the  number 
of  coupled  states  appearing  in  (1.6).  In  a process 
involving  multiple  surfaces  we  usually  need  to  include 
those  states  that  are  "close"  to  the  initial  surface,  in 
the  sense  that  they  are  energetically  accessible;  we  often 
need  to  include  closed  channels  (energetically 
inaccessible  states)  as  well.  Hence  we  are  often 
confronted  with  the  need  for  excited- state  surfaces;  these 
are  usually  not  of  the  same  quality  as  the  lower  lying 
states. 

An  alternative  to  this  approach  is  to  consider  as  our 
electronic  basis  states  which  are  not  eigenstates  of  H 

el 

This  notion  dates  back  as  far  as  Zener  (1932),  and  the 
term  diaba tic  to  describe  a state  which  does  not 
adiabatically  adjust  to  the  instantaneous  position  of  the 
nuclei  was  coined  in  comparatively  recent  times  (lichten, 
1963;  1 967)  . It  is  interesting  to  note  that  until  this 
time,  the  prevalent  notion  cf  molecular  theorists  was  that 
the  diagonal  states  were  the  electronic  states. 
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Onf ortunately , there  is  no  unique  set  of  states  which 
satisfy  the  requirement  of  a diatatic  representation;  they 
lack  qeneral  definition.  The  later  work  of  Smith  (1969) 
was  significant  in  this  connection.  Smith's  representation 
consisted  of  choosing  the  electronic  basis  in  such  a 
manner  that  the  off-diagonal  terms  T"  and  T*  vanish;  in  a 
sense,  the  Born-Oppenheimer  approximation  becomes  exact  in 
this  representation.  Variations  on  this  theme  have  these 
terms  small  but  not  necessarily  zero.  The  reader  is 
directed  to  the  recent  work  of  Delos  and  Thorson  (1979)  in 
this  connection.  In  all  such  pictures,  however,  the 
dominant  coupling  appears  in  the  V terms,  the 

ij 

off-diagonal  matrix  elements  of  H . The  large  gradients 

el 

found  in  the  adiabatic  wa vef unct ions  become  small  or,  in 
the  case  of  Smith's  representation,  zero.  Diabatic  states 
are  not  constrained  to  obey  a non-crossing  rule  and  go 
smoothly  through  "avoided  crossings."  The  non-crossing 
rule  is  thus  not  a law  of  nature  but  appears  only  as  an 
artifact  of  a particular  choice  of  representation  of  the 
electronic  states  of  molecules. 

This  completes  our  survey  of  the  mathematical 
background.  We  will  now  review  a few  of  the  theoretical 
approaches  to  collisional  ionization  in  the  next  section. 
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1.3  Survey  of  Theories  of  Collisional  Ionization 
Ionizing  collisions  in  atoms  of  the  type 
A*  + B — >A+B++e~ 

— > AB+  + e- 

vhere  A is  electronically  excited  have  been  the  subject 
of  considerable  experimental  and  theoretical  interest  as 
these  processes  may  have  large  reactive  cross  sections 
even  at  low  (thermal)  energies.  For  the  case  where  the 
electronic  excitation  energy  of  A exceeds  the  ionization 
potential  of  B,  the  initial  discrete  N-electron  electronic 
state  is  degenerate  with  the  continuum  of  electronic 
states  corresponding  to  the  remaining  N-1  bound  electrons 
and  the  free  electron.  Schematic  potential  energy  curves 
for  N and  N-1  particle  states  are  displayed  in  Figure  1; 
the  lower  curve  is  the  threshhold  of  an  electronic 
continuum.  The  initial  (A*,B)  state  may  thus  be  thought 
of  as  an  autoicnizing  or  resonance  electronic  state. 
Among  the  important  theoretical  aspects  of  this  process 
are  the  coupling  between  discrete  and  continuum  electronic 
states,  the  coupling  among  continuum  electronic  states, 
and  finally  the  influence  of  these  coupled  states  upon  the 
heavy  particle  dynamics. 

Classical  treatments  (Miller  and  Schaefer,  1970)  are 
based  on  a close  encounter  of  metastable  and  target,  with 
autoionizaticn  occurring  at  the  classical  turning  point 
(distance  of  closest  approach)  in  the  effective  (AB*  plus 
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centrifugal)  potential.  Further,  the  kinetic  energy  of 
relative  motion  of  A*  and  B must  exceed  the  maximum  in  the 
effective  potential  for  reaction  to  occur.  A slightly 
more  detailed  description  {Miller,  1970)  makes  use  of  an 
autoionization  width  Y for  decay  of  the  discrete  state 
into  the  continuum.  Transitions  may  occur  at  separations 
greater  than  the  classical  turning  point,  so  a transition 
probability  may  be  obtained  for  the  heavy  particles  both 
on  the  inward  and  outward  parts  of  their  trajectory. 

Semiclassically  one  obtains  a probability  amplitude 
(or  S-matrix  element)  as  a superposition  of  probabilities; 
differential  and  total  cross  sections  then  depend  upon  the 
sguare  modulus  of  the  amplitude  and  so  differ  from  the 
classical  result  in  that  interference  terms  (arising  from 
interference  between  transitions  during  in  and  out 
trajectories,  for  example)  are  new  present. 

Since  the  resonance  state  lies  in  the  continuum  of 
the  ionized  channel,  quantum  mechanical  treatments  have 
been  based  largely  upon  a complex,  local  adiabatic 
potential  which  governs  the  heavy  particle  dynamics 
(Nakamura,  1969;  Miller,  1970;  Miller  and  Schaefer,  1970). 
This  complex  potential  results  in  complex  phase  shifts, 
and  the  imaginary  part  of  the  potential  accounts  for  loss 
of  flux  (via  inelastic  and/or  reactive  processes)  from  the 
initial  channel  (Mott  and  Massey,  1971;  Sicha,  1974).  A 
formal  development  involving  an  effective  Hamiltonian 
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derived  using  projection  operator  algebra,  and  which 
contains  the  local,  complex  potential  as  an  approximation, 
will  be  presented  in  detail  in  Chapter  II.  However,  at 
this  point  we  wish  to  review  a recent  and  novel 
theoretical  approach  (Bellum  and  flicha,  1978)  to 

chemiionization  to  set  the  stage  for  both  the  developments 
in  Chapter  II  and  a portion  of  the  numerical  results  to  be 
presented  in  a subsequent  chapter. 

For  the  particular  collisional  ionization  processes 

we  wish  to  study.  Penning  and  associative  ionization,  the 

final  ionized  channel  consists  of  bound  cr  unbound  nuclei 

plus  a free  electron.  Consequently,  an  electronic  basis 

consisting  solely  of  discrete  (bound)  electronic  functions 

is  inappropriate  here.  Since  the  outgoing  electron  can 

have  arbitrary  kinetic  energy  (constrained  only  by 

conservation  of  total  energy  during  the  collision),  we 

need  to  augment  the  discrete  basis  by  a continuum  of 

electronic  states  <J>  (6)  ; here  6 labels  the  asymptotic 

kinetic  energy  of  the  free  electron.  A prescription  for 

constructing  such  resonance  and  continuum  electronic 

states  is  given  in  the  next  chapter.  Assuming  for  the 

moment  that  only  a single  initial  discrete  state  <J>  is 

d 

significant,  we  expand  the  total  wavefunction  for  nuclei 

and  electrons  (in  a body-fixed,  center  of  mass  of  the 

, BF-CMN)  in  the  basis  [$  , { (6)  } ] : 

d 


nuclei  frame 
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f = <J>  X + S (0,  inf)  d€  <J>  {€)  X (€)  , 

d d 


(1.16) 


Implicit  in  the  above  is  a neglect  of  nuclear  and 

electronic  angular  momentum  coupling  and  the 

representation  of  the  outgoing  electron  by  a pure  s-wave 

(a  more  detailed  discussion  will  be  presented  in  Chapter 

II) . In  what  fellows  we  shall  restrict  out  attention  to 

PI;  some  remarks  concerning  extension  of  the  theory  to  AI 

will  be  made  at  the  conclusion  of  this  section. 

Substitution  into  the  SchrOdinger  equation,  taking  the 

usual  scalar  products  and  making  a partial  wave  expansion 

for  the  nuclear  motion  leads  to  a continuous  set  of 

coupled,  second-order  ordinary  differential  equations  for 

the  heavy  particle  functions  X , X (€)  . 

d 

To  make  the  problem  computationally  tractable,  Bellum 

and  iiicha  (1978)  introduced  a complete,  orthonormal  basis 

of  discrete  functions  of  the  continuum  variable  e, 

ff  (6)}.  Expanding  the  continuum  coefficients  in  this  set 
J 


X (€)  = 2(J=1,N)  f (€)  X 

J J (1.17) 

and  proceeding  in  the  usual  manner  gives  rise  to  a matrix 
SchrOdinger  equation  of  order  N+1  describing  motion  of  the 
heavy  particles  with  angular  momentum  partial  wave  L : 

L L L L 

D X = U X (1.18) 


where 


17 


L 

(D  ) 

ij 


(d2/dr2+k  2-L  (L+1) /R2)S 
i ij 


(1.19) 


for  i , j = 1 to  N+ 1 , 


L 

(U  ) = (2M/X2)V 

1 1 d 


L 

(0  ) = (2 M/)<2)  (V  + e ) 

j j + J 

for  j=J+1=2  to  N + 1,  and 

L 1/2 

(U  ) = (2M/X2)  {/J  e)  V 

i j d J 


(1.20) 

(1.21) 


(1.22) 

for  i=1,j=J+1=2  to  N+1  and  j=1,i=J+1=2  to  N+1,  and  zero 

otherwise.  Channel  wavenumbers  and  the  reduced  mass  are 

given  by  k and  M.  The  potentials  V and  V are  diabatic 
i d + 

potential  energies  for  the  initial  discrete  state  and 

molecular  ion  state,  respectively.  The  coupling  V is  the 

dJ 

electronic  matrix  element  between  the  discrete  state  and 

the  J-th  discretized  continuum  state;  e is  the  kinetic 

J 

energy  of  the  ejected  electron  in  the  J-th  interval  of  the 
continuum.  The  choice  of  the  set  is  that  of  Bellum  and 
flicha  (1978),  and  is  repeated  here;  for  J = 1,2,...,N  : 

-1/2 


f (6) 

J 


iA  6) 


(1.23) 
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for 


e - (4  €)  /2  < € < € ♦ (4  6)/2 

J J 


and  zero  otherwise.  An  important  feature  to  note  is  that 

the  couplings  between  continuum  electronic  states 

(V  ,J*K)  have  been  neglected.  The  potential  matrix  thus 
JK 


blocks  as  follows: 


V 

V 

• • • 

V 

d 

Id 

2d 

rd 

d 1 

V +6 
+ 1 

0 

• • • 

0 

d2 

0 

v +e 
+ 2 

• « • 

0 

V 

Nd 

0 

0 


V 0 0 

dr 


V 0 0 

dN 


V +6 
♦ r 

0 


The  boundary  condition  for  (1.18)  at 
0)  reads 


0 


(1.24) 


V +6 

+ N 
_ 

the  origin  (R  = 


L 

X (0)  = 0 (1.25) 


and  from  the  asymptotic  behavior  of  the  wavefunction  we 
obtain  the  real  and  symmetric  reaction  matrix  K (or  B 
matrix)  : 


L L L L 
lim  X = J ♦ N K 

R — > oo 


(1.26) 


where  J and  N are  diagonal  matrices  of  Biccati  Bessel 
functions  of  the  first  and  second  kind  (Abramowitz  and 


Stegun,  1972) 
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L c -1/2 

(J  ) = a k J (k  F) 

ij  ij  j L j (1.27) 

L c -1/2 

(N  ) = o k N (k  P)  . 

ij  ij  i L j (1.28) 

Recall  that  we  have  restricted  our  treatment  to  PI,  so 

only  open  channels  (k  2 > 0)  of  nuclear  motion  are 

1 

considered.  The  unitary  and  symmetric  scattering  matrix 

S is  given  in  terms  of  K (Mott  and  Massey,  1971)  by 

L L L 

S = (1  + iK  ) (1  - iK  )-»  (1.29) 

The  obvious  advantage  to  using  K matrix  boundary 

conditions  is  that  we  may  work  entirely  with  real 

quantities,  rather  than  the  complex  functions  dictated  by 

S matrix  boundary  conditions.  Within  the  above  treatment, 

the  partial  wave  reactive  scattering  amplitude  is 

L L 

S (E)  = 2.(J=  1,  N)  f (£)  S (E) 

Cl  J J1  (1.30) 

L 

where  the  S are  amplitudes  for  collision  with  emission 
J1 

of  an  electron  with  kinetic  energy  distribution  given  by 

the  probability  amplitude  f (6) . The  Penning  ionization 

J 


cross  section  per  unit  energy  of  the  emitted  electron 
(also  known  as  Penning  Electron  Energy  Distribution  or 


PEED)  is  thus 
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-1 

d(T  (E,6)  /d6  = n/k  2 G) 

d 

L 

*2(L)  ( 2L  ♦ 1 ) |S  (E)|2 

J1  (1.31) 


and  the  Penning  ionization  cross  section  is  given  by 

6"  (E)  =5(0,e')  d€  d6"(E,e)/d6  (1.32) 

where  €'  is  the  maximum  value  o£  the  e mi tted-electron 
kinetic  energy  for  PI  and  is  fixed  by  energy  conservation 
once  the  initial  collision  energy  is  specified. 

In  the  previous  work  (Bellum  and  Micha,  1978)  the 
neglect  of  the  continuum-continuum  couplings  suggested  the 
following  approximation:  replace  the  N+1  - dimensional 

matrix  problem  by  N 2-state  problems  (in  essence,  each 
discrete-continuum  transition  is  taken  to  be  independent 
of  all  ethers) , namely 


L L L L 

D Y = U Y 

J J J J (1.33) 

L 

where  we  explicitly  write  ()<2/2M)U  as 

J 


1/2 


V 

d 


(A  6)  V 

dJ 


(1.34) 


1/2 

(A  6)  v v + e 

dJ  + J 

for  J=1  to  N.  We  will  examine  in  detail  (by  way  of 
numerical  studies)  the  validity  of  this  decoupling 
approximation  in  a later  chapter,  and  in  particular  will 
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focus 

on 

the  cros 

s sections 

(1.31) 

and 

(1. 

32) 

resulting 

from 

the 

full 

matrix  pr 

oblem 

and 

from 

the 

N 2-state 

problems. 

In 

view  of 

the  an 

ticipated 

reduction  in 

computational  effort  in  solving  N pairs  of  coupled  second 

order,  ordinary  differential  equations  versus  a matrix 

differential  eguation  of  order  N+1,  it  is  critical  to  have 

knowledge  of  the  range  of  applicability  of  this 

approximation;  this  has  implications  for  both  this  and 

other  coupled  channels  treatments  of  PI  (including  recent 

work  on  PI  in  a laser  field:  Bellum  and  George,  1978; 

1979;  Bellum,  Lam  and  George,  1978)  ). 

When  we  wish  to  extend  the  coupled  channels  theory  to 

the  case  of  associative  ionization,  where  the  heavy 

particles  remain  bound  after  electron  emission,  we 

discover  at  once  that  the  boundary  conditions  of  this 

problem  are  such  as  to  complicate  the  solution.  For  closed 

channels,  k 2 < 0 and  the  matrix  elements  of  J and  N are 

1 

given  in  terms  of  modified  spherical  Bessel  functions  of 

the  first  and  third  kinds  (Abramowitz  and  Stegun,  1972) 

L r 1/2 

(J  ) =6  (k  R)  I (k  R) 

ij  ij  j L+1/2  j (1.35) 

L c 1/2 

(N  ) =5  (k  R)  K (k  R)  . 

ij  ij  j L*1/2  j (1.36) 


The  K matrix  then  blocks  as  follows  : 
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K K 

oo  oc 

K K 

CO  cc 


(1.37) 


where  K , K , K and  K are  the  open-open, 
oo  oc  co  cc 

open-closed,  closed-open  and  closed-closed  sub-matrices  of 

K.  The  S matrix  is  then  expressed  in  terms  of  K in 

oo 

analogy  with  (1.29).  However,  the  modified  Bessel 

function  K pertaining  to  closed  channels  exhibits 

L+  1/2 

exponentially  decreasing  behavior  with  increasing  R,  and 
we  realize  that  there  will  be  numerical  difficulties  in 
trying  to  extract  the  amplitude  of  this  decaying 

component. 

Apart  from  this  obstacle  which  precludes  extension  of 
this  coupled  channels  methodology  to  AI,  there  are  other 
problems  of  a fundamental  nature  when  one  attempts  to 
apply  the  discretization  concept  in  a general  case.  One 
difficulty  is  to  reconcile  the  general  notion  of  a 
discretization  of  the  electronic  continuum  with  the  fact 
that  the  final  nuclear  states  in  AI  are  discrete  (i.e., 
the  vibration-rotation  levels  of  a molecular  ion).  In 
other  words,  a portion  of  a discretization  basis  ought  to 
span  the  discrete  part  of  the  electron  kinetic  energy 
range  (corresponding  to  AI)  as  well  as  the  continuum  part 
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{corresponding  to  PI) . This  suggests  that  we  need  to  know 
in  advance  the  detailed  vibrational-rotational  structure 
of  the  molecular  ion.  In  PI,  the  same  discretization 
basis  was  used  for  each  heavy  particle  angular  momentum 
partial  wave.  By  contrast,  for  AI  the  basis  would  need  to 
be  different  for  each  partial  wave:  a different  set  of 
vibrational  levels  (with  different  level  spacings)  is 
accessible  for  each  L value.  Furthermore,  if  we  wished  to 
generalize  the  preceding  model  to  include  higher  angular 
momentum  partial  waves  for  the  emitted  electron,  this  may 
introduce  yet  another  angular  momentum  dependence  of  a 
discretization  basis  set.  One  last  aspect  of  this  problem 
is  the  functional  form  of  the  discretization  basis;  the 
set  of  functions  of  Bellum  and  Micha  (1978)  and  used  by 
Relyea  and  Micha  (1979)  was  chosen  on  physically  appealing 
grounds.  For  a more  general  choice  of  basis  (e.g., 
overlapping  functions)  the  potential  matrix  U does  not 
block  as  in  (1.24)  and  the  possibility  of  invoking  an 
analogous  (and  expedient)  decoupling  approximation  may  be 
lost. 

Although  the  numerical  study  suggested  by  the  above 
remarks  is  beyond  the  scope  and  intent  of  the  present 
work,  a final  assessment  of  this  theory  must  await 
numerical  experiments  as  regards  discretization  basis  set 
type  and  size,  as  well  as  an  examination  of  the  validity 
of  other  approximations  that  have  been  made.  The 
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computational  scheme  foe  the  solution  of  the  FI  part  of 
the  collisional  ionization  problem  will  be  presented  in 
Chapter  III,  together  with  detailed  results  pertinent  to 
the  Penning  ionization  of  Ar  by  He*(1s2s,3S)  at  thermal 
energies. 

In  the  next  chapter  we  will  present  an  alternative 
approach  to  the  coupled  channels  treatment  described  in 
this  section.  It  is  an  effective  Hamiltonian  theory  based 
upon  projection  operator  techniques,  and  we  will  see  that 
many  of  the  deficiencies  of  the  coupled  channels  method 
are  remedied  in  the  new  theory.  In  a later  chapter  we 
propose  methods  of  solution  of  the  resulting  scattering 
equations  at  various  levels  of  approximation. 
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CHAPTER  II 

EFFECTIVE  HAMILTONIAN  FORMALISM  FOR 
IONIZING  ATOMIC  COLLISIONS 

II.  1 Electronic  Wavefu ncti ons  Appropriate  to  Penning__and 

Associative  Ionization 

Before  proceeding  to  develop  a theory  for  the  heavy 
particle  dynamics  cf  collisional  ionization  in  atomic 
systems,  consideration  must  first  be  given  to  the 
electronic  wavef unctions  appropriate  to  such  processes. 
Once  a prescription  has  been  given  for  the  construction  of 
electronic  wa vef unctions  for  the  initial  AB*  and  final 
(AB+, e~)  states,  potential  energy  functions  and  couplings 
may  then  be  used  in  a theory  of  heavy  particle  dynamics. 

As  has  been  already  mentioned  in  the  preceding 
chapter,  the  initial  discrete  electronic  state  of  AB*  is 
embedded  in  a continuum  of  states  of  the  type  (AB*,e~)  and 
as  such  is  an  autoionizing  electronic  state,  capable  of 
decay  by  electron  emission  even  in  the  limit  of  zero 
relative  velocity  of  the  nuclei.  Interpreted  as  a 

molecular  autoionization  for  fixed  nuclei,  the  transition 
rate  would  be  independent  of  the  B-0  couplings  which 
dynamically  couple  electronic  and  nuclear  motions.  This 
transition  will  be  governed  by  a coupling,  through  an 
electronic  Hamiltonian,  of  the  (as  yet  unspecified) 
discrete  and  continuum  electronic  states;  for  the  low 
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energy  collisions  to  be  considered  in  this  work  it  is 
reasonable  to  assume  that  such  couplings  still  dominate 
the  process.  These  remarks  suggest  that  we  must  construct 
a diabatic  ( ncn-diago na 1)  representation  of  the  electronic 
Hamiltonian;  we  will  now  review  some  conventional 
technigues  for  obtaining  such  wavefunctions  with 
particular  attention  to  the  description  of  the  continuum 
electronic  state. 

Following  the  historical  development  where  the 

application  first  was  to  three-electron  systems  (He*  * H) , 

Miller  et  al.  (1972)  chose  to  represent  both  the  initial 

discrete  function  $ and  the  molecular  ion  function  <t>  by 

i + 

configuration  interaction  (Cl)  wavefunctions.  For  ♦ this 

i 

is  an  application  of  the  "stabilization"  method  of  Hazi 
and  Taylor  (1970).  It  has  long  been  argued  that  although 
autoicnizing  or  resonance  states  are  not  true  electronic 
eigenfunctions  (in  that  they  require  some  constraint  for 
their  definition),  they  can  still  be  obtained  by  the 
established  procedures  for  calculating  ordinary  bound 
states:  one  lcoks  for  a positive  energy  root  which  is 
stable  with  respect  to  basis  set  change  and  identifies 
such  a state  as  a resonance  cr  "bound  state  in  the 
continuum."  A continuum  electronic  state  is  constructed  by 
augmenting  the  determinants  in  by  a continuum  orbital 
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of  energy  6;  this  form  ensures  a properly  antisymmetrized 
continuum  electronic  state.  For  the  continuum  orbital  a 
Coulomb  orbital  centered  on  the  hydrogen  atom  and 
orthogonalized  to  all  the  bound  MO's  was  employed.  A 


Coulomb  orbital  is 

a 

reasonable 

description  for 

an 

electron 

outside 

of 

an  ionic 

core;  the  effect 

of 

distortion 

at  shorter 

distances 

is  neglected. 

The 

autoionization  width 

is 

given  by 

K= 

2nj>|<<>  (€ ; 1, 

2,. 

.,N)|H  -E | 

el 

♦ (1,2,-.,N)>|2 
i (2. 

D 

where  p is  the  density  of  final  states  and  H the 
J el 

electronic  Hamiltonian.  The  wavefunction  <J>  (6  ; 1 , 2, . . , N)  is 

a continuum  function  in  which  the  free  electron  has  energy 

6 and  orientation  € in  the  tody-fixed  frame.  It  was  in 

anticipation  of  the  width  calculation  that  the  same  set  of 

MO's  was  used  for  both  c)>  and  $ ; this  eliminated 

i + 

the  need  for  taking  matrix  elements  between  determinants 
containing  non-orthogonal  orbitals.  Bound-free  matrix 
elements  were  evaluated  numerically.  Technical  details 
aside,  this  approach  ensures  the  orthogonality  of  bound 
and  continuum  states  and  guarantees  that  the  continuum 
state  obeys  the  boundary  condition  appropriate  to  a 
scattering  (positive  energy)  state  at  large  distances. 

A more  recent  formulation  (Hickman  et  al.,  1977), 
for  He*  ♦ H,  represented  the  molecular  ion  state  $ 


again 
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by  a single  determinant  and  constructed  N-electron  states 
by 

* (1,2,...,N)  = A { $ (1,2,...,M-1)U  (N)} 
m + m (2.  2) 

where  (u  } is  some  MO  basis  and  A is  the  antisymmetrizer. 
m 

These  new  states  have  the  physical  form  of  an  ionic  core 

plus  another  orbital;  they  are  not  continuum  states. 

Defining  a projection  operator  P by 

P = Zltn=1  ,M)  |1>  {1,2,...  ,N)  ><«>  (1,2,...,N)  | 

m m (2.3) 

they  diagonalized  (1-P)H  ( 1 — P)  in  some  multi-determinantal 

el 

basis  (Cl)  to  extract  those  bound  states  orthogonal  to 
states  in  the  P subspace.  The  ground  state  of  HeH  was  not 
one  of  the  eigenvalues  of  QHQ  and  was  obtained  by 

diagonalizing  the  full  electronic  Hamiltonian.  All  states 
of  QHQ  were  thus  autoionizing  or  resonance  states  of 
positive  energy;  the  desired  eigenvalue  was  the  lowest 


root  of  QHQ. 

The  resonance 

state 

energies 

were 

essentially 

identical 

with  those 

obtained  from 

the 

stabilization 

technique 

applied 

to 

the  full 

H , 
el 

demonstrating 

the  advantage  to 

diagona lizing 

only  a 

sub-block  if 

a prudent 

choice 

of 

P is  made. 

It 

also 

appears  that  one  key  to  the  success  of  this  method  is  to 
include  enough  configurations  in  P so  that  the  resonance 
is  the  lowest  solution  of  QHQ  but  yet  to  leave  Q "large 


enough"  to  give  an  accurate  resonance  function.  A 

continuum  wavefunction  <)>  (6 ; 1, 2,  . . . , N)  was  constructed  as 

4>  (€;  1 ,2,  . . . ,N)  = Afd)  (1,2,...,N-1)4»(6;U)} 

♦ (2-4) 

where  <}>  (0 ; N>  is  a Coulomb  orbital  orthogonalized  to  the 

orbitals.  An  important  aspect  of  this  work  was  the  neglect 

of  the  non-square-integrable  portion  of  the  orthogonalized 

Coulomb  orbital.  Justification  for  this  procedure  was 

found  in  the  fact  that  for  a set  of  MO's  which  is  a good 

basis  in  a localized  region  in  space,  the  continuum  part 

of  this  orthogonalized  continuum  orbital  will  be  large 

only  outside  this  region.  Since  the  dominant  contribution 

to  the  auto  ionization  width  came  from  the  inner  region, 

this  neglect  did  not  constitute  a severe  numerical  error. 

Although  computationally  expedient,  the  "continuum"  state 

thus  obtained  now  obeys  boundary  conditions  appropriate  to 

bound-state  wa vef unctions,  namely  asymptotic'  decay. 

Further,  the  resonance  and  continuum  electronic  states 

were  orthogonal  only  because  the  same  MO  basis  was  used  in 

the  construction  of  <j>  , and  <J>{6;N). 

i + 


Another  scheme  for  obtaining 

autoionizing 

or 

resonance 

states 

of  molecules  is 

the  calculation 

of 

so-called 

Siegert 

eigenvalues  (Siegert,  1939);  in 

more 

recent  applications  (Bardsley  and  Junker,  1972;  Isaacson, 
1979;  Isaacson  and  Miller,  1979;  McCurdy  and  Rescigno, 
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1979)  orbital  basis  sets  for  Cl  are  augmented  by  complex 
orbitals  asymptotically  representing  Coulomb  orbitals,  and 
the  resulting  complex  Hamiltonian  matrix  is  diagonalized 
to  obtain  eigenvalues  and  vectors.  Applications  of  any  of 
the  above  techniques,  capable  of  producing  very  high 
quality  electronic  wavef unctions  and  ionization  widths, 
are  in  practice  limited  to  few-electron  systems;  treatment 
of  larger  systems  at  the  Hartree-Fock  or  post-HF  level  is 
ordinarily  facilitated  by  using  atomic  core 
pseudopotentials,  restricting  excitations  when 
constructing  Cl  determinants,  etc. 

For  the  ionization  of  Ar  by  He*(1s2s,3S)  we  will  need 
an  electronic  wavefunction  for  (or  the  potential 
corresponding  to)  the  initial  3Z.  state;  for  the  final 
channel  molecular  ion  state  which  separates  to  He  + 
Ar*(2P)  we  would  in  principle  need  both  2Z.and  2n  states. 
As  our  ultimate  interest  will  be  in  collision  energies  on 
the  order  of  0.  1 eV,  we  may  safely  neglect  additional 
states  in  either  channel.  The  closest-lying  state  of 
higher  energy  for  HeAr*  corresponds  to  He(1s2p,3P)  * Ac 
and  lies  (asymptotically)  over  1 eV  (Moore,  1949)  above 
the  initial  state;  an  excited  state  of  HeAr*  corresponding 
to  He  + Ar* (3s 1 3p6, 2S)  lies  over  13  eV  above  He  ♦ 
Ar*  (3s23p5,  2P)  . 

Assuming  the  projection  of  electronic  orbital  angular 
momentum  is  conserved  in  the  collision  process,  we 
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conclude  that  the  emitted  electrons  are  in  either  C-/  or  n 
orbitals.  However,  recent  experimental  evidence  (Hoffman 
and  Morgner,  1979)  lends  strong  support  to  the  idea  that 
only  <T  electrons  are  emitted;  i.e.,  HeAr*  is  formed  only 
in  the  2Z  state.  This  may  at  first  be  surprising  in  that 
the  statistical  weights  for  forming  2Z  and  2n  are  in  the 
ratio  of  1:2  in  collisions  of  He  and  Ar*.  However,  we  must 
bear  in  mind  that  in  PI  and  AI  the  molecular  ion  is  formed 
not  by  collision  of  He  and  Ar*  but  by  decay  of  HeAr*  by 
electron  emission.  On  the  basis  of  the  preceding  remarks 
we  are  justified  in  treating  the  heavy  particle  dynamics 
as  a two-state  process:  a 3Zlstate  for  the  initial  channel 
and  a 2Z1  state  for  the  molecular  ion  channel. 

For  the  initial  "hole"  state  ( 3Z)  , inspection  of 
separ ated-ato m energies  of  He,  He*,  Ar  and  Ar*  (Moore, 
1949)  indicates  that  it  is  not  the  lowest  energy  state  of 
that  symmetry;  for  the  open-shell  ion,  HeAr*  (2Z)  is  the 
lowest  state  of  its  symmetry.  For  recent  results  on  HeAr* 
the  reader  is  directed  to  the  work  (SCF,  Cl)  of  Olson  and 
Liu  (1978)  and  to  the  SCF  results  of  Sidis  et  al.  (1977). 
Although  no  ab-initio  results  are  available  for  HeAr*,  a 
one-electron  model  has  recently  been  developed  by  Siska 
(1979a).  fiesults  for  both  systems  at  the  MSX^  level  have 
been  previously  reported  by  Bellum  and  Micha  (1974;  1977). 
Our  effort  does  not  focus  just  on  obtaining  initial  or 
molecular-ion  wa vefunctions  of  arbitrary  accuracy  but 
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rather  is  directed  at  using  such  results  to  obtain 
wavef unctions  and  potentials  appropriate  to  collisional 
ionization  processes.  Our  requ irements  on  the 
wavef unctions  are  1)  that  the  discrete  and  continuum 
states  are  orthogonal  (recall  out  close-coupling  expansion 
in  Chapter  I was  in  terms  of  an  orthogonal  basis)  and  2) 
that  we  have  a realistic  description  of  the  continuum 
orbital,  since  in  low-energy  collisional  ionization  "slow" 
electrons  are  ejected,  and  these  will  have  interaction 
with  the  ionic  core  which  deviates  from  Coulombic  at  short 
distances. 

The  major  difficulty  in  obtaining  a proper 
description  of  a continuum  state  is  in  the  representation 
of  the  outgoing  electron.  For  very  high  energy  electrons 
a plane  wave  (PW)  or  orthogonalized  plane  wave  (OPW) 
treatment  is  cften  adequate.  At  lower  energies  Coulomb 
wave  (C W)  or  orthogonalized  Coulomb  wave  (OCW)  are  better 
descriptions  as  they  represent  the  true  field  experienced 
by  the  electron  at  large  distances  from  the  ion.  These 
approaches  all  share  the  same  defect:  an  assumption  for 
the  form  of  the  free  electron  wavef unction.  Since  this 
orbital  represents  a scattering  state  it  is  perhaps  more 
natural  to  approach  its  description  using  the  tools  of 
scattering  theory,  although  fcr  collisional  ionization 
this  approach  seems  to  have  been  overlooked.  Simply 
stated,  we  could  attempt  to  solve  the  problem  of  the 
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scattering  of  a low-energy  electron  by  an  open  shell 
molecular  ion.  Indeed,  in  the  fixed-nuclei  approximation 
a formulation  has  been  given  for  electron  scattering  by 
closed  shell,  neutral  targets  by  Burke  and  Sinfailam 
(1970)  (in  the  static  exchange  approximation)  and  by 
Chandra  and  Tetnkin  ( 1976)  (for  a static  plus  model 
polarization  potential) . Thus  the  tools  are  available  for 
the  construction  of  reliable  potentials  with  which  to 
solve  the  electron-molecule  (and  maybe  molecular  ion) 
scattering  problem.  However,  this  constitutes  a rather 
ambitious  undertaking  and  we  would  prefer  to  utilize 
traditional  bound  state  techniques  to  improve  upon  and 
extend  the  work  discussed  at  the  beginning  of  this 
section. 

We  now  briefly  outline  a scheme  whereby  the 

wavef unctions  (whose  matrix  elements  are  the  desired 
potentials)  may  be  constructed. 

1}  First  obtain  a molecular  ion  wavefunction 

<*  (1,2,.. .,8-1)  for  the  open  shell  (2X)  state  which 


+ 

asymptotically  correlates  with  He  + Ar*(2P).  As  it  is  the 
lowest  of  Z symmetry  it  is  likely  to  be  obtainable  by 
variational  means,  although  a post-HF  method  may  be  needed 
to  ensure  correct  dissociation. 

2)  To  go  beyond  the  Coulomb  approximation  for  the 
emitted  electron  a straightforward  approach  is  proposed 
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here.  We  choose  to  adopt  the  method  of  Harris  (1967) 
which  has  been  applied  extensively  to  atomic  problems. 
First  we  obtain  a wavefunction  for  the  molecular  ion  HeAr+ 
at  the  HF  level.  In  general  the  solution  will  yield  both 
bound  (negative  energy)  and  victual  (unoccupied)  molecular 
orbitals.  The  virtual  orbitals  are  not  stationary  (except 
in  special  circumstances)  with  respect  to  basis  set 
variation  as  they  do  not  correspond  to  bound  electronic 
states.  These  orbitals  can  be  thought  of  as  providing  a 
representation  for  the  scattering  of  an  electron  by  HeAr+ 
in  its  22Z  state.  A continuum  election  orbital  can  be 
expanded  in  a L2  basis  augmented  by  regular  and  irregular 
non-L2  functions  which  properly  account  for  scattering 
boundary  conditions;  a purely  L2  expansion  can  only 
provide  a good  description  of  a scattering  function  in  a 
localized  or  truncated  region  cf  space  and  does  not 
satisfy  correct  boundary  conditions.  For  an  electron 
outside  a neutral  atom  or  molecule  these  scattering 
functions  could  be  regular  and  irregular  spherical  Bessel 
functions;  for  a molecular  ion  the  proper  choice  would  be 
Coulomb  functions.  Application  of  the  Harris  method 
proceeds  as  follows:  at  incident  electron  kinetic  energies 
equal  to  each  of  the  positive  virtual  orbital  energies 
and  for  the  choice  of  the  HF  orbitals  as  the  L2  basis,  the 
ratio  of  the  amplitudes  of  the  non-L2  functions  is 
obtained,  and  this  defines  the  scattering  orbital 
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wavef unction.  Modifying  the  computational  basis  yields  a 
new  set  of  victual  orbitals  with  new  energies,  and 
repeating  the  above  procedure  we  can  generate  a collection 
of  continuum  orbitals  at  practically  any  energy.  This 
treatment  is  often  termed  a static-exchange  treatment; 
strictly  speaking,  exchange  is  really  only  included  in  the 
target  wavef unction.  To  explicitly  include  exchange 

effects  between  the  incident  electron  and  the  target 
requires  antisymmetciza tion  of  the  total  electronic 
wavef unction ; this  effect  can  be  less  rigorously  included 
by  addition  of  an  approximate  exchange  potential  in  lieu 
of  antisymmetrization.  Effects  of  polarization  of  the 
molecular  ion  by  the  free  electron  can  also  be  included  by 
model  polarization  potentials.  Approximations  to  exchange 
and  polarization  contributions  are  common  techniques  in 
current  electron-molecule  coupled  channel  calculations. 
We  can  now  construct  a collection  of  continuum  electronic 
wavef unctions  of  the  form 

<M6;  = [♦  (1f2,...,N-1)#(€;H)  ] 

+ 


where  <HS;N)  is  the  scattered  orbital;  € represents  the 
orientation  6 and  magnitude  € of  the  emitted  electron 
energy  in  the  body-fixed  frame. 

3)  To  obtain  the  molecular  resonance  state  we  could 
first  calculate  a HF  ground  state  { 1 , 2 , . . . , N)  for  HeAr 
(i^).  A zeroth  order  approximation  to  the  hole  state 
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( 1,  2,. , N)  could  be  obtained  by  a particular  single 
hs 

excitation  out  of  This  wavefunction  is  then 

orthogonalized  to  the  continuum  w avef unctions  obtained  in 

2)  to  give  a crude  description  of  the  molecular  resonance, 

$ . An  improved  treatment  could  be  obtained  by 
r 

constructing  a projector  P in  the  continuum  function  space 
of  2)  and  solving  the  Cl  problem  for  (1-P)  H ( 1-P)  - The 
basis  could  be  the  orbitals  of 

hs 

The  above  scheme  yields  wavefunctions  which  satisfy 

the  requirements  of  a diafcatic  representation;  since 

neither  resonance  nor  continuum  states  are  eigenfunctions 

of  the  full  electronic  Hamiltonian,  they  will  have 

non-zero  matrix  elements  with  H . It  also  satisfies  our 

el 

two  previously  stated  requirements  on  the  continuum 
wavef unction : orthogonality  to  the  resonance  state  and  an 

improved  description  of  the  outgoing  electron.  It  is 
hoped  that  the  remarks  here  may  provide  stimulation  and 
motivation  to  determine  quantities  of  interest  in 
scattering  theory.  For  diatomic  systems  of  reasonably 
liqht  atoms  where  detailed  dynamical  studies  are  possible, 
accurate  potentials  and  couplings  ace  desirable  since  then 
an  accurate  assessment  of  the  scattering  theory  may  be 


made. 
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II. 2 Effective  Hamiltonian  Formalism 
With  electronic  states  obtained  (in  principle)  using 
the  methods  of  the  previous  section,  we  will  now  formulate 
the  heavy  particle  scattering  problem  to  which  these 
states  are  needed  input.  Based  upon  the  remarks  in 
Section  1.3  such  a reformulation  is  necessary  if  the 
theory  is  to  be  able  to  describe  both  Penning  and 
associative  ionization  processes.  The  methods  of 
Feshbach's  unified  theory  of  nuclear  reactions  (1958, 
1962)  will  be  applied  to  the  problem  of  ionizing 
collisions;  the  method  of  projection  operators  allows  the 
many-channal  (infinite  in  number  in  the  present  work) 
problem  to  be  recast  in  terms  of  an  effective  Hamiltonian 
for  one  of  the  channels  of  interest. 

The  input  to  our  treatment  consists  of  a set  of 

Born-Oppenheimer  electronic  states;  one  of  these  is  a 

discrete  state  $ which  corresponds  to  the  molecular 

d 


resonance  state  of  the  previous  section,  and  which 
correlates  asymptotically  with  separated-atom  limits  of  A* 
f B-  In  light  of  the  anticipated  application  we  will  have 
a single  continuum  of  electronic  states  {<)>  (6)  } associated 
with  a particular  state  of  the  molecular  ion  AB*  and  which 
correlates  asymptotically  with  the  limits  of  A ♦ B+  + e~; 
here  6 denotes  both  magnitude  e and  orientation  e of  the 
electron  in  the  body-fixed  frame.  Additional 


AB  + 
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electronic  states  give  rise  to  multiple  continua,  each 
having  a different  lower  limit  or  threshhold;  the  theory 
is  generalizable  to  this  situation  but  the  development 
here  will  focus  on  the  one  discrete  - one  continuum  case. 
The  only  restrictions  upon  the  electronic  basis 

r (cj>  (6) } ] are  the  orthogonality  and  normalization 

d 

conditions 

<<|>  | <t>  > = 1 

d d (2.5) 

<4>  li>  (e)  > = o 


d 

(2.6) 

V 

= S (IT'-e)  . 

(2.7) 

<<M6')  I4>(6)> 

Since 

in 

the  present 

work  we  are 

interested  in 

e luci 

datin 

g the 

effects  of 

nuclear  motion 

on  electron 

emiss 

ion. 

the 

viewpoint  is 

fundamentally  d 

ifferent  from 

that 

of  Na 

kamura 

(1959;  1971) 

and  Fujii  et 

al.  (1970), 

where  Fano's  method  (1961)  for  autoionizing  atomic  states 
was  extended  to  molecular  problems.  In  these  works  a 
totally  adiabatic  framework  was  chosen;  configuration 
interaction  of  discrete  and  continuum  electronic  states 
resulted  in  a level  shift  and  width  for  the  resonance 
energy.  This  energy  is  necessarily  a local  quantity  since 
the  theory  was  based  upon  a fixed-nuclei  description  of 
the  ionization.  This  complex  resonance  energy  was  then 
used  as  the  potential  for  nuclear  motion;  decay  of  the 
initial  state  by  electron  emission  was  governed  by  the 
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local,  imaginary  part  of  the  potential.  The  theory  yielded 
easily  obtainable  ionization  cross  sections  but  is 
incapable  of  providing  insights  into  the  interplay  of 
electronic  and  nuclear  motion  which  is  of  interest  here. 

Our  method  is  to  construct  an  operator  which  governs 
heavy  particle  motion  in  the  incident  (A*  + B)  channel  and 
which  incorporates  effects  due  to  both  direct  (elastic) 
and  reactive  (ionizing)  collisions.  To  accomplish  this  a 
projection  operator  Q is  defined  onto  the  molecular 
resonance  state: 


Q = |(f  ><4>  | . 

d d (2.8) 

As  defined  here  Q projects  only  onto  the  space  of 
electronic  wa vef unctions  and  has  as  a consequence  only  a 
parametric  dependence  upon  the  icternuclear  coordinate  B. 
Our  choice  of  Q is  similar  to  that  of  O’Malley  (1966)  in  a 
study  of  dissociative  attachment,  where  a projector  was 
chosen  onto  the  closed  electronic  channels.  The  Q 
channel,  although  closed  with  respect  to  electronic 
motion,  is  an  open  channel  for  the  nuclei.  The  operator 
Q,  together  with  its  complement  P 


P = 1 - Q (2.9) 

are  the  basis  for  the  effective  Hamiltonian  formalism  to 
be  presented  here.  To  the  extent  that  other  electronic 


states  can  be  neglected,  P becomes  simply 

P = 5 (0  , inf  ) d€  | (6)  ><$  (6)  | 


(2. 10) 
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which  implies  nothing  more  than  that  the  electronic  basis 

constitutes  a complete  set  for  our  purposes.  The 

SchrOdinger  equation  partitions  as 

0 (E-H)  Qf  = QHPf 

Q Q (2.11) 

P ( E-H)  Pf  = PHQf 

Q Q (2.12) 

where  we  have  used  Q as  a subindex  on  the  wavefunction  to 
denote  that  the  initial  channel  is  in  the  Q subspace;  here 
H is  the  full  Hamiltonian.  We  also  introduce  solutions  to 
the  corresponding  homogeneous  equations 

Q (E-H)  Qf  = 0 

q (2.13) 

P (E-H)  Pf  = 0 . 

P (2.14) 

With  the  initial  channel  as  Q,  the  formal  solution  is 
given  by 

P|«-=G  + PHCt»  + 

Q P Q (2.15) 

Qf  *=QV ♦ + G + QHPV  + 

Q q Q Q 

= Qf  + + G+QHPG+PHQf ♦ 

q Q p Q 

=Qf+  + G+F*Qf+ 

q Q Q (2.  16) 

where 

F*=QHPG*PHQ 

P (2.17) 

and  the  reduced  Green’s  functions  are 


G+=[  Q (E ♦ - 

Q 


H)Q  ]-» 


(2. 18) 
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and 


G*=[P(E+  - H)  P ]~ 1 . 

P (2.19) 

The  effective  Hamiltonian  for  initial  channel  motion  is 

obtained  by  noting  that  for  H given  by  (1.2)  and  for  the 

choice  of  Q,  the  Q-space  projection  of  the  total 

scattering  wavefunction  for  nuclei  and  electrons  is 

Ql*=  I4>  ><<*  IH*> 

Q d d Q 

= 14  >|X*> 

d Q (2.20) 

and  then 

(E  - V - T - <<t>  | F+ 1 <)>  >)  J X4’>=0 

d 8 d d Q (2.21) 

where 


V =<0  | H | g>  > . 

d d el  d (2.22) 

The  physical  content  of  (2.21)  is  the  following:  the 

interaction  of  A*  with  B is  given  by  the  potential  V ; 

d 

this  potential  accounts  only  for  elastic  scattering  in 

this  channel  and  does  not  allow  for  inelastic  or  reactive 

processes.  The  term  <<J>  IF+I4  > is  responsible  for 

d d 

ionization  through  the  coupling 

PHQ  = J d6|  <>(€)  ><<M€)  |H  | <t>  ><♦  | 

el  d d (2.23) 

followed  by  propagation  in  the  ionized  (AB*,e~)  channel 

governed  by  G*  and  then  possible  transition  through  QHP 
P 
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back  to  the  initial  A*  + B state.  The  heavy  particle 
motion  in  the  initial  channel  is  thus  governed  by  an 
'•optical"  potential  V 


V = V + <4.  |F+M>  > 

d d d (2.24) 

whose  components  give  rise  to  direct  and  reactive 

scattering.  it  is  clear  from  (2.17)  and  (2.19)  that  the 

potential  is  complex  and  energy  dependent.  For  a 

discussion  of  effective  Hamiltonians  for  molecular 

collisions  the  reader  is  directed  to  the  work  of  Micha 

(1969;  1974;  1976)  and  the  review  of  Rabitz  (1976).  From 

this  point  forward  our  treatment  departs  from  that  of 

Feshbach  (1958;  1962)  and  it  is  worthwhile  to  point  out 

the  distinction  in  this  work.  Following  Feshbach  we 

rewrite  the  second  term  in  V as 


<$  |QHP[P(E+  - H)P]-tPHQ)4>  > 

d d (2.25) 

and  note  that  the  magnitude  of  this  term  will  vary  rapidly 

when  E is  in  the  neighborhood  of  an  eigenvalue  of  PHP. 

Assuming  for  the  moment  a set  of  discrete  eigenstates  {X  } 

i 

of  PHP,  with  eigenvalues  (E  } , we  write  (4>=<J)  ) 

i d 

V=V  + Z(i*j)  CcMQHPI  X >(E  + -E  )~1<X  |PHQ|<J>> 
d i i i 

+ < <J>  I QHP  | X >(S+-E  ) ~ 1 <X  |PHQ|<|>> 

j 3 j (2.26) 

where  E is  an  eigenvalue  of  PHP  in  the  vicinity  of  E.  For 
3 
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a given  E we  have  written  V as  a sum  of  three 

contributions,  the  first  two  of  which  presumably  vary 

slowly  with  E and  which  give  rise  to  a smooth  E-dependence 

of  the  cross  section;  for  E neat  E we  would  anticipate  a 

j 

large  contribution  to  the  cross  section,  a phenomenon 
which  is  loosely  given  the  term  "resonance."  The  foregoing 
remarks  have  been  applied  to  a single,  isolated  resonance 
and  must  be  modified  (Feshbach,  1958;  1962;  Lemmer,  1966) 
when  there  are  several  eigenvalues  (isolated  or 

overlapping)  of  PHP  in  the  vicinity  of  E. 

For  our  application  the  splitting  of  the  optical 
potential  V does  not  offer  any  practical  advantage.  The 
initial  state  in  PI  and  AI  is  itself  a resonance  (of  a 
different  kind)  . That  the  initial  channel  energy  may 
coincide  with  an  eigenvalue  of  PHP  is  of  secondary 
importance  since  resonances  of  the  type  illustrated  in  the 
Feshbach  treatment  have  net  been  observed  in  PI  or  AI  (a 
different  kind  of  resonance  behavior  has  been  uncovered  in 
PEED  using  the  model  of  Section  1.3,  and  further 
discussion  of  such  will  be  deferred  until  Chapter  III). 

For  our  choice  (2.10)  of  P we  obtain 

g + = 3 de  3d6*  i<MG)>g+(€,e*)  «>(€•)  I 

P (2.27) 

and 

F+  = Jde  5 d€*  V _,g  + (€,€•) 

d€  ed 


(2.28) 
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where 


and 


g*(e,e')  = <<M6)  | (E+-fl  -T  }-t  | ♦(€•)> 

el  R 


(2.29) 


V = <4>  |H  |(M€)  > . 

d el  (2.30) 

The  operator  g*  is  an  operator  in  the  space  of  heavy 

particle  functions.  At  this  point  we  introduce  two 

approximations  which  together  will  considerably  simplify 

the  operators  G*  and  F*  : assume  that  1)  we  have  an 

P 


approximately  diagonal  representation  of  H within  the  P 

el 

subspace,  and  2)  the  continuum  electronic  energy  is  that 

of  the  N-1  particle  state  plus  the  asymptotic  energy  6 of 

the  outgoing  electron;  we  have 

<♦(€•)  I H |<M6)>  = (V  €)S  (€•-€) 

el  ♦ (2.31) 

where 


V =<4  (1,2,.  .,N-1)  |H  (1,2,..,  N-1)  10  (1,2,..  ,N-1)  > 

* + el  ♦ (2.32) 

and  g a molecular-ion  wavef unction.  This  approximation 
+ 

was  employed  by  Bellum  and  Micha  (1978)  to  render  the 
matrix  second-order  differential  equation  more 
computationally  tractable;  recall  that  the  decoupling 
approximation  (Section  1.3)  immediately  followed  this 
step.  Me  note  here  that  such  further  approximation 
(decoupling)  will 


not 


be  needed  in  the  present 
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formulation;  this  will  become  evident  in  Chapter  IV. 

Physically,  the  neglect  of  continuum-continuum  couplings 

(2.31)  stems  from  the  fact  that  the  ejected  electron 

leaves  the  region  of  the  molecule  rapidly  compared  to  the 

heavy  particle  relative  velocity;  there  is  little 

likelihood  of  transitions  between  continuum  states 

corresponding  tc  different  electron  kinetic  energies. 

With  this  we  have  the  diagonal  forms 

G*  = 5 36  | <t>  {€)  >g  + (€)  <<J>  (6)  | 

P,d  (2.33) 

and 

F + = j d€  V _.g+  (6)  V_^ 
d de  Cd  (2.34) 

and  where 

g+  (6)  = (E*  - T - (V  «•  6)}~i  . 

R ♦ (2.35) 

We  thus  obtain  an  inhomogeneous  differential  equation 

for  the  Q projection  of  the  total  scattering  wavefunction 

{E-V  (R ) — T JX+ftT,^)  = 
d R Q 

3dR'  jd6  V _^(R)  <R|  g«-  {€)  |R'>V^  (R ' ) X*  (E,  it*  ) 

d€  €d  Q (2.36) 

where  we  have  assumed  spherical  symmetry  of  the  initial 

channel  potential  V and  of  the  coupling  V _J>.  To  proceed 

d d£ 

with  the  analysis  we  will  formally  expand  g*  (6)  in  its 
spectrum  of  bound  and  continuum  heavy  particle  states.  As 


we  may  ultimately  wish  to  focus  separately  on  either 
Penning  or  associative  ionization  we  will  treat  each  of 
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these  cases  in  turn. 


For  PI  introduce  a complete  set  of  heavy  particle 
functions  in  the  ionized  channel 

1 = J I * ("E')  ><r  (?')  I 

PI  (2.37) 


( E • > 0 and  continuous)  where 

(E'-T  -V  ) | 1 (E')  > = 0 

R ♦ (2.38) 

and  where  1 denotes  the  unit  operator  in  this  subspace 
PI 

of  the  total  heavy  particle  function  space.  Making  a 
partial  wave  expansion 


K'E'/B)  =2(1,  in)  Y (E')Y  (R)i1 

lm  lm 

♦exp  (iS  ) V ( E ' , R)  /R 

1 1 (2.39) 


we  obtain 

<B  | q*  (€)  | S 1 > = 1/(RR*)£(l#n)  ) dE'  Y (R) 

PI  lm 


♦ Y (R*)  ? (E',B)f  (E  ' , R* ) 

lm  1 1 


*(E*  - 6 - E ' } ~ 1 (2.40) 

where  the  functions  i|  (E',R)  satisfy  the  radial  equation 

1 

{d2/dR2+ (2M/#2)  (E*-V  (R) ) -1  (1+1) /R2}W  (E',R) 

+ 1 


= 0 (2.41) 

subject  to  the  asymptotic  boundary  condition 

lim  1 (E',R)  = A (/)  sin  (k'R-ln/2  + S ) 

B — >oo  1 1 


(2.  42) 
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where  S is  the  real  partial  wave  phase  shift  and  A(k)  = 
1 

X-V2n-^2. 

For  AI  introduce  the  complete  set  of  bound  heavy 
particle  states 

N n N n 

1 = Z(nJN)  | V(E  ) ><V(E  ) | 

AI  J J (2.43) 

n 

(E  <0  and  discrete)  where 
J 

n N n 

(E  -T  -V  ) |f  (E  ) > * 0 . 

J B + J (2.44) 

Expanding  in  partial  waves 

N n _i  /v  n 

|(E  , R)  = Y <B)  ?{E  ,B)/B 

J JN  J J (2.45) 

we  obtain 

<B|g+(6)|P‘>  a 1/ (RP'  )£(n,J,N)  Y (E) 

AI  JN 

a * n n * 

*Y  ( R ' ) ? (E  ,B)f  (E  ,8') 

JN  J J J J 


n 

*(E*  - 6 - E }-» 

J (2.46) 

n 

where  the  functions  }}  (E  ,B)  are  solutions  to  the 

J J 

vibration-rotation  eigenvalue  equation 

n 

(d2/dH2  - (2M/K?)  V (E)  - J (J+  1 ) /E2  ) V <B  , R) 

♦ J J 
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n n 

= (2H/*T2)  E V (E  ,F)  . 

J J J (2.47) 

flaking  a partial  wave  expansion  of  the  initial 

channel  heavy  particle  function  X + 

Q 

-s.  -*■  * A .A 

X * (E,R ) =21.(1.,  M)  Y (E)Y  (R) 

C Lfl  LM 


*i  exp  (id  ) X + (E,R)  /R 


L L 

(2.40) 

substitution  into 

(2.36),  using  (2.40) 

and  (2.46), 

A 

A * 

pre  mu  It i plica  tion 

by  (2M/)(2)  Y (E)  Y 

(R)  and 

L*  S ’ L ' 

M* 

A 

A 

integrating  over  dE 

and  dR  gives 

L 

(d*/dR2  + (2M/K2)  ( E~ V (R)  -L  (L+1)  /R*}  i 

d 

*exp  <iS  ) X (E,  R)  * 

L L 


(2H/K2)  JaE^dR  Y (E)Y  (R)  $ d R ' J d€  V ^(R) 

LN  LN  d€ 


*{  I + II  }V^  (S')  I(L',}I')  Y (E)  Y (R) 
ed  L ' N r L ' N1 


*i  exp  (iS  ) X (E,B')/R' 

L'  L‘  (2.49) 

where  the  quantities  I and  II  are  given  by  (2.40)  and 

(2.46),  respectively.  Asymptotically, 

lim  X (E,R)  = A (k)  sin  (kR-ln/2  + S ) 

B-->oo  l L (2.50) 
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where  c5  is  the  complex  partial  wave  phase  shift. 
_L 

Alternatively  for 

S = + i§2 

L L L (2.51) 

X = X»  + i X2 

L L L (2.52) 

we  may  write 

X » = A(k)coshS2  sin(kF  - Ln/2  * &l) 

L L L (2.53) 

X2  = A(k)sinhS2  cos(kR  - Ln/2  + . 

L L L (2.54) 


The  phase  shift  is  complex  because  the  potential  is 
complex,  the  imaginary  part  of  which  gives  rise  to  flux 
loss  through  ionization. 

To  further  simplify  the  derivation  we  need  to  examine 

the  form  of  the  coupling  matrix  element  V of  (2.30).  We 

dS 


assume  that  the  continuum  orbital  4>{6;N)  in  the  continuum 

wavef unction  $ (6;  1,2,  • . . N)  can  be  expanded  in  angular 

momentum  partial  waves  (in  the  body-fixed  frame)  as 

<Me;N)  = <f  (e,  r ) 

N 


= Z.  { 1,  m ) Y (€)  <#  (€lm  ;r  ) 

lm  N (2.55) 

and  we  need  to  determine  how  many  partial  waves  1 should 

be  retained  in  the  expansion.  To  the  extent  that 

rotational  Dorn-Oppenheimer  couplings  can  be  neglected, 

the  most  straightforward  approach  is  to  proceed  as  in 
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Hickman  et  a 
compute  part 
and  deterrain 
have  apprec 
resonance 
projection 
angular  mome 
conservation 
projection . 


1.  (1977)  and  Miller  et  al.  (19  72)  and  simply 
ial  widths  at  various  internuclear  separations 
e the  minimal  number  of  partial  waves  1 which 
iable  contribution.  If  the  symmetries  of  the 
and  molecular  ion  states  ace  known,  the 
m of  the  emitted  electron  electronic  orbital 
ntum  on  the  internuclear  axis  can  be  fixed  by 
of  total  electronic  ornital  angular  momentum 
That  is. 


ra  = M*  - M* 

where  M*  (M*j  is  the  projection  for  th 
(molecular  ion)  state.  The  tacit  assumption 
the  axial  component  of  the  electronic  or 
momentum  is  a good  quantum  number.  In  the 
Micha  and  Nakamura  (1975)  analyzed  molec 
correlation  diagrams  in  the  separated-  and 
limits  for  the  initial  and  ion  states  to  ded 
number  of  partial  waves  1 which  should  be  inc 
for  HeAr*,  HeAr*  this  set  was  1=0, 1,2 
dominant  contribution) . This  work  was  exten 
1976;  Bellura  and  Micha,  1977)  to  includ 
Born-Op penheimer  rotational  couplings 
contributing  partial  waves  for  several  a 
atom-ion  pairs.  Since  heavy  particle  angular 
large  (typically  cn  the  order  of  10  to  100) 
electron  angular  momentum,  we  are  justified  i 


e resonance 
here  is  that 
bital  angular 
same  spirit 
ular  orbital 
united-atom 
uce  a minimal 
luded  (e. g. , 
with  1=0  the 
ded  (Bellum, 
e effects  of 
to  deduce 
tom-atom  and 
momentum  is 
compared  with 
n neglecting 
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electronic  angular  momentum  in  both  the  initial  and  final 
channels.  As  it  is  the  object  of  the  present  work  to 
focus  on  a detailed  description  of  the  heavy  particle 
dynamics  in  ionizing  collisions,  we  can  with  justification 
simplify  our  task  by  considering  only  the  first  (1=0)  term 
in  the  expansion. 

Our  previous  result  (2.49)  becomes 

{d2/dR2+  (2M/X2)  (E-V  (R) ) -l  (L  + 1)/R2)  X (E,R}  = 

d L 

(2M/Jf2)  jdR'  J d€  V (R)[  JdE'  |f  (£‘,R) 

€d  L 

* 

*f  (E , R ' ) {E^-G-E*} “* 

L 

n n * 

+ 2.(n)  f (E  , R)  V (E  ,S») 

L L L L 


n 

*{E+-€-E  }-»]  V (R  * ) X (E  ,R  * ) 

L d€  L (2.56) 

and  where  the  right-hand  side  becomes,  on  evaluating  the 

quantity  inside  [...], 

R.H.S.  = (2M/X2)  j dR  •[  5 dE'  f (E’,F>  {p.v.JdG 


V (R)  V (R  * ) (E-G  -E')-»  - in  V (R) 
d€  €d  dE  * ' 


*v  (R*)}  V (E 1 , R ' ) 
E»  ’ d L 


♦ Z(n) 


f (E  ,3 ) {p.  V.  ) de  V (P) 

l l de 
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n 

*V  (R ' ) (E-e  -E  )-i  - in  V (P) 

€d  L dE'" 

n * 

*V  (RM  } I (E  , R ' ) ] X (E,RM 

E'"d  L L L (2.57) 

n 

vfhere  E''=E-E'  and  E,,,=E-E  and  p.v.  denotes  a Cauchy 

L 

principal- value  integral. 

The  first  two  terras  in  [...]  in  (2.57)  are 

contributions  to  the  potential  from  continuum  (Penning) 

channels  and  the  last  two  terms  are  from  AI  channels.  The 

result  (2.57)  clearly  demonstrates  the  non-locality  and 

energy  dependence  of  the  optical  potential.  Also  note 

that  all  discrete  and  continuum  final-channel  heavy 

particle  states  of  a given  L are  coupled  to  the  initial 

channel  motion  of  the  same  angular  momentum  partial  wave 

L.  We  symbolically  write  the  last  result 

0 X (E,R)  = (2M /k2)  jdfl'  K (E;R,R')X  (E,R  ') 

L L L L (2.58) 

where 

D = (d*/dR2+ (2M/XP)  (E-V  (P)  ) - L (L+  1 ) /R 2 } 

L d (2.59) 

and  the  kernel  K is  decomposed  as 

I 

PI  AI 

K (E  ; R r R ' ) = K ( E ; R , R ' ) + K (E  ; R f R 1 ) 

L L L (2.60) 

or 


K ( E ; R , R ' ) = A (E ;R ,R  ' ) 
L L 


(i/2)  Y (E;  R, R ' ) 
L 


(2.61) 
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and  where 

PI  AI 

A { E ; R , R 1 ) - A (E  ; R , R • ) + A (E; R, R' ) 

L L L (2.62) 

and 


PI  AI 

Y (E;RrR»)=K  (E;R,R')+y  (E ; R,  R ' ) . 

L E L (2.63) 

In  contrast  to  the  set  of  coupled  second-order  ordinary 

differential  equations  we  obtained  in  Section  1.3,  in  the 

present  theory  we  have  a single  inhomogeneous  differential 

equation  for  an  "effective"  nuclear  wavef unction ; the 

complexity  of  the  problem  lies  in  the  kernel  K (E;R,R ' ) . 

L 

To  complete  the  formalism  we  need  the  scattering 

cross  sections  for  PI  and  AI.  Cross  sections  for 

transitions  from  states  i to  f ate  given  as  the  transition 

rate  R per  unit  incident  flux;  R is  given  by  (Micha, 
fi  fi 

1974) 

R = (2tt/X)  |T  |2^(E  -E  ) 
fi  fi  f i (2.64) 

where  T is  the  relevant  transition  amplitude.  For 
f i 

continuum  states  normalized  to  unity  on  the  energy  scale 

<»(E*)  lf(E)>  * S(E»-E)  (2.65) 

the  flux  is  kVC  (2n ) *X]- 
i 

For  Penning  ionization  the  final  channel  variables 
are  E'  and  €;  from  the  results  of  Appendix  B we  have 
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T (E,E',8)  = <|I-(E*)V  X*(5)> 

fi  ed  Q (2.66) 

where  <...>  denotes  JdR.  Expanding  and  X in  partial 

waves  ( (2.39)  and  (2.48),  respectively)  we  obtain  the 

triply  differential  cross  section 

i \ /" 

d^E,  E' ,6) /dE'de  = cd(E  -E  )|T  |2 

i f fi  (2.67) 

where 

c = (2n/K)  [ (2n)  3j</k2  ] 

i (2.68) 

and 

IT  |2  = |Z(LM)  Y (E')Y  (E)exp(iS  ) 
fi  LM  LM  L 

*exp  (io  )jdR  |-(E',R)V  (F)  X+  (E , R)  | 2 . 

L L Gd  L (2.69) 

The  energy  distribution  of  emitted  electrons  is 

obtained  by  integrating  over  heavy  particle  angles  and 

energies  and,  by  choosing  the  initial  relative  momentum  to 

lie  along  the  z-axis  in  the  center  of  mass  system,  we 

obtai n 


PI 

d(T  (E,G)/dG  = c|Z(L)  [ ( 2L+  1 ) / ( 4n)  ] V2 

*exp  (iS  ) exp  (ic)  ) 

L L 

*jdR  R)V  (R)  X*  (E,R)  | 2 . 

L €d  L 

The  Penning  ionization  cross  section  is  simply 

PI  f PI 

^ =j  (0,e»)  de  d6“(E,6)/de 

with  8*  the  maximum  value  of  the  electron 


(2.70) 


(2.71) 
for  PI, 


energy 
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determined  by  energy  conservation. 

For  associative  ionization  the  final  channel 


variables  are  E and  € and  we  have  the  amplitude 
J 

n n _> 

T (E,E  ,6)  =< | ( E ) V X ♦ (E)  > . 
fi  J J €d  Q (2.72) 

Making  partial  wave  expansions  of  the  heavy  particle 

functions  we  arrive  at 

c n 

IT  I2  = [ ( 2 J+ 1 ) /4tt  ] |)dR  f (E  ,P.) 
fi  * J J 

*V  (R)  X+  ( E , R)  exp  (iS  ) | * 

ed  j j 


and  the  differential  cross  section  per  unit  energy  is 

AI  n £ 

d6"  (E,E  ,€)/d€  = co  (E  -E  )|T  |*  . 

J i f fi  (2.73) 

The  cross  section  for  formation  of  the  molecular  ion  in 

the  vibrational-rotational  state  (n,J)  is  obtained  by 

integrating  (2.73)  over  €,  which  gives 

AI  n n 

6“  (£»E  ) =c[  (2J+ 1) /4n  ]|)  dR  \ (E  ,R) 

J J J 


*V  (R)  X+  ( E,  R)  exp  (iS  ) | 2 

G"d  J j (2.74) 

where  in  the  last  result  we  have  made  use  of  the 

energy-conservation  relation 


E+^=€"+E 


with  ^ the  asymptotic  energy  separation  between  resonance 
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and  molecular  ion  potential  curves.  The  AI  counterpart  to 

the  electron  energy  distribution  for  PI  is  a discrete  set 

of  lines  for  the  electron  energy  spectrum.  For  a 

particular  vibrational  state  n we  have 

AI  Ain 

(E,n)  = Z(J)  <T  (ErE  ,6) 

J (2.75) 

and  the  total  AI  cross  section  (for  formation  of  all 
possible  molecular  ion  vibrational  and  rotational  states) 
is 

AI  AI  n 

C (E)  = 2(nJ)  C (E,  E ,S) 

J (2.76) 

This  completes  our  theoretical  development  of  an 
effective  Hamiltonian  formalism  for  ionizing  atomic 
collisions.  In  the  next  chapter  we  will  discuss  some 
computational  aspects  of  the  coupled  channels  theory  of 
Section  1.3  and  present  detailed  numerical  studies  of  PI 
of  Ar  by  He(1s2sf3S).  In  the  final  chapter  we  will  return 
to  the  effective  Hamiltonian  theory  where  we  present  a 
method  of  solution  for  the  inhomogeneous  second-order 
ordinary  differential  equation  (2.58). 


CHAPTER  III 

COMPUTATIONAL  ASPECTS  OF  A COUPLED  CHANNELS  THEORY  FOR  PI 
III.  1 Solution  of  the  Matrix  Schrodinqer  Equation 
la  order  to  obtain  numerical  solutions  to  a set  of  N 
coupled  second-order  ordinary  differential  equations 
(ODE’s)  such  as  (1.18)  we  must  first  specify  2N  constants. 
Regularity  of  solutions  at  the  origin  (1.25)  serves  to 
determine  N of  these  constants,  and  the  remaining  N can  be 
chosen  as  the  values  of  the  first  derivatives  at  the 
origin.  The  specification  of  the  N-dimensional  derivative 
matrix  is  otherwise  arbitrary  so  long  as  it  leads  to 
linearly  independent  solutions.  Suitable  choices  are  the 
unit  matrix,  a diagonal  matrix  whose  elements  have  ’’small" 
values  (such  that  solutions  start  out  with  snail  slopes) 
or,  more  generally,  any  non-singular  matrix. 

Choosing  the  origin  as  a starting  point  is  frequently 
a source  of  difficulty,  however.  Regions  interior  to 
classical  turning  points  of  interatomic  potentials  are 
non-classical  regions,  and  it  is  in  these  regions  that  the 
solutions  decay  exponentially  with  decreasing  R.  As  a 
consequence,  in  outward  integration  from  the  origin  one  is 
often  plagued  by  an  exponential  growth  of  seme  of  the 
solution  components  and  an  attendant  loss  of  linear 
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i ndep  enJence  of  the  v ect  or  s.  Two  reined  ies  suggest 
themselves:  modify  the  boundary  conditions  at  the  origin 
or  periodically  orthogona lize  the  solution  matrix.  The 
first  choice  could  be  effected  by  re- specif yi ng  the  origin 
at  some  value  of  R intermediate  between  the  origin  and  the 
classical  turning  point;  the  choice  of  new  origin  could  be 
by  trial-and-er cor  (Bernstein,  1960)  or,  following  Bellum 
and  Micha  (1978)  , by  using  semiclassical  approximations  to 
the  wavefunction  and  derivative  and  selecting  as  the 
origin  that  value  of  R at  which  the  slope  is  less  than 
some  specified  (usually  small)  tolerance.  In  a 
multichannel  problem  one  may  encounter  a situation  in 
which  the  classical  turning  points  are  widely  separated, 
so  that  the  origins  for  each  channel  will  likewise  be 
separated.  During  the  early  stages  of  outward 
propagation,  however,  the  exponential  growth  problem  may 
still  be  present  unless  one  suppresses  those  solutions 
with  the  outermost  origins.  The  presence  of  centrifugal 
barriers  (also  nonclassical  regions)  in  any  channel  may 
again  tend  to  introduce  linear  dependencies  requiring 
periodic  orthogonalization  of  the  solutions  in  these 
areas . 

For  the  Penning  ionization  problem,  a computer 
program  DEVOG  was  written  to  solve  a set  of  coupled 
seccnd-order  ODE's.  The  algorithm  of  deVogelaere  (1955), 
generalized  by  Lester  (1968,  1971)  to  the  multichannel 
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(matrix)  case,  was  selected  as  the  integrator.  This  is  a 
fourth-order  method  as  is  comparable  in  accuracy  to 
fourth-order  Runge-Kutta  (RK)  methods.  As  an 

error -controlled  mesh  selector  (where  step-length  is 

modified  to  satisfy  a local  error  criterion)  was  not 
developed  for  the  deVogelaere  method  until  the  recent  work 
of  Coleman  and  Mohamed  (1979) , a fixed  step— length 
algorithm  was  used.  Step  size  was  chosen  to  be  dependent 
upon  the  collision  energy;  typically  10  steps  per 
deBroglie  half- wavelength  were  found  to  be  adequate.  A 
single  channel  (elastic  scattering)  option  of  the  program 
was  satisfactorily  tested  against  scattering  phase  shift 
results  of  Bernstein  (1960)  who  employed  a standard 
R unge-K u tta- Gil 1 (RKG)  integrator;  the  coupled  channels 
option  was  successfully  tested  against  partial  wave  S 
matrix  elements  of  Olson  and  Smith  (1971;  1972),  where  a 
Numerov  integration  scheme  was  used. 

For  application  to  Penning  collisions  we  have  made  a 

study  of  the  ionization  of  Ar  by  He*(1s2s,3S)  at  thermal 

energies.  A two-state  approximation  has  been  invoked:  the 

He*  and  Ar  are  assumed  to  move  in  a 32L  potential,  and  the 

motion  of  the  heavy  particles  He  and  Ar*  is  governed  by  a 

2Z.  potential;  these  are  the  potentials  V and  V , 

d * 

respectively,  appearing  in  (1.34).  The  functional  forms 
of  these  potentials,  as  well  as  the  coupling  matrix 
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elements  V (1.34)  were  those  used  in  previous  work 
dJ 

(Bely ea  and  Micha,  1979;  Bellum  and  Micha,  1978))  and  are 
given  in  Appendix  A.  It  should  te  pointed  out  that  while 
"realistic"  functional  forms  are  chosen  for  the 

potentials,  the  origins  of  the  parameters  lie  both  in 
theory  and  experiment  and  so  to  this  extent  serve  as 
models  with  which  we  will  do  detailed  dynamical  studies. 
In  this  connection,  the  reader  is  directed  to  the  recent 
work  of  Siska  (1979a)  and  Clson  and  Liu  (1978)  for  a 
discussion  of  32l  and  2Z.  potentials,  respectively.  At  a 
collision  energy  of  65  meV  (0.0024  Hartrees)  we  have 
computed  the  PEED  (1.31)  within  the  decoupling 

approximation  (1.33)  of  Chapter  I.  The  step  size  was 
chosen  to  be  1/20-th  of  the  smallest  deBroglie  wavelength 
in  each  of  the  coupled  channels.  The  discretization  basis 
step  size  was  £ 6 = 2.0(10-*)  a.u.  for  6 < 0.  1503888 

a.u.,  and  2J  e = 0.5(1 0—  ♦ ) a.u.  for  6 greater  than  this 
value.  In  Figure  3 the  results  are  presented  for  a 
maximum  L value  of  50  and  a discretization  basis  of  35 
functions.  The  maximum  L value  was  selected  on  the  basis 
cf  convergence  of  the  Penning  ionization  cross  section 
upon  inclusion  of  higher  partial  waves.  The  continuum  mesh 
was  decreased  for  the  upper  range  of  electron  kinetic 
energy  to  more  clearly  elucidate  the  structure  suggested 
in  the  PEED  in  earlier  work  (Bellum  and  Micha,  1 978).  In 
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all  cases  the  symmetry  of  the  K matrix  (1.26)  was  better 
than  1 part  in  103.  Total  computation  time  to  solve  the 
ccupled  equations  (1.33)  and  calculate  cross  sections  was 
54  minutes  on  an  IBM  370/165  computer;  this  is  roughly  2 
seconds  per  partial  wave  and  per  6,  and  was  faster  than 
the  computations  of  the  earlier  work  by  a factor  of 
one- third. 

The  qualitative  features  (e.g.,  overall  shape  and 
position  of  the  maximum)  are  in  general  agreement  with  the 
experimental  results  of  Hotop  (1974).  That  there  is  more 
structure  in  the  computed  results  is  not  surprising;  the 
widths  of  the  peaks  in  the  computed  spectrum  are  on  the 
order  of  10~s  a.u.  ( < 1 meV) , whereas  experimental 
spectra,  limited  by  thermally  broadened  beams  and  finite 
detector  resolution,  have  a full-width-half-maximum 
resolution  of  10  - 20  meV  (Hotop  and  Niehaus,  1969).  One 
is  led  to  question  whether  any  interpretation  can  be  given 
to  the  structure  of  the  calculated  distribution  or  if  the 
maxima  are  merely  computational  artifacts. 


As 

we  saw  in  the 
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E + & = € + E ; 

f (3.1) 


here  E and  6 are  measured  relative  to  the  separated-atom 
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energies  of  the  colliding  species  in  the  entrance  channel, 

and  ^ and  E are  measured  relative  to  the  final  channel 
f 

molecular  ion  asymptote.  The  interpretation  of  the  single 
feature  in  the  PEED  of  Bellum  and  Micha  (1978)  was  based 
upon  the  notion  of  an  orbiting  resonance  in  the  final 
channel;  the  colliding  particles  may  become  temporarily 
bound  by  a centrifugal  barrier,  increasing  the  likelihood 
of  reaction  during  this  prolonged  encounter.  The  final 
channel  energy  was  speculated  to  lie  just  above  a maximum 
in  the  centrifugal  barrier  for  L = 24,  the  partial  wave 
which  dominated  the  cross  section  at  that  peak.  In 
addition,  the  classical  turning  point  in  the  initial 
channel  was  found  to  lie  directly  above  this  barrier.  The 
combination  of  these  two  effects,  the  increased 
probability  of  a transition  in  the  vicinity  of  a classical 
turning  point  and  the  possibility  of  some  sort  of  resonant 
state  in  the  final  channel,  was  thought  to  explain  the 
enhancement  in  the  spectrum  at  that  particular  electron 
energy. 

We  have  pursued  this  line  of  questioning  by  studying 
partial  wave  contributions  to  the  PEED  at  each  value  of 
the  emitted  electron  kinetic  energy  where  there  is  cross 
section  enhancement;  these  are  labelled  (a)  - (g)  in 
Figure  3.  One  first  observation  was  that  the  classical 
turning  point  for  the  dominant  partial  wave  in  the 
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entrance  channel  was  located  above  (but  not  necessarily  at 
the  iraxiinum  of)  a centrifugal  barrier  for  that  same 
L-value  in  the  final  channel,  consistent  with  the 
observation  in  the  earlier  work.  To  facilitate  further 
discussion.  Table  1 contains  dominant  L-value 
contributions  to  the  PEED,  as  well  as  final  channel 
energies  and  partial  wave  barrier  maxima.  For  all 
structures  except  (e)--the  largest  one-- the  final  channel 
energies  lie  just  above  the  barrier  maxima;  these  may 
correspond  to  scattering  resonances,  where  the  heavy 
particles  in  the  final  channel  have  just  barely  enough 
energy  of  relative  motion  to  avoid  binding.  These  states 
may  be  loosely  thought  of  as  "bound  states  in  the 
continuum"  of  nuclear  motion,  and  this  opens  an 

interesting  area  for  further  study  utilizing  techniques  of 
molecular  electronic  structure  theory. 

The  so-called  "stabilization"  method  of  Hazi  and 
Taylor  (1970)  is  based  upon  the  observation  that  in 
basis-set  electronic  structure  calculations  one  obtains, 
from  an  N-dimensicnal  matrix  eigenvalue  problem,  perhaps  M 
bound  (negative  energy)  states  and  N-H  positive  energy 
states.  By  adjusting  the  basis  the  bound  levels  will 
remain  stationary  (assuming  the  basis  is  "good  enough"  to 
begin  with)  while  the  energy  of  the  scattering  states 
fluctuates;  in  the  event  that  there  are  "virtual"  levels 
whose  energies  do  not  change  with  basis  set  adjustment. 
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these  levels  are  said  to  have  been  "stabilized"  and  are 
termed  resonances,  bound  states  in  the  continuum  or 
autoionizing  electronic  states.  Turning  to  the  nuclear 
motion  problem,  the  same  basis  set  techniques  can  be 
applied  to  determine  the  vibrational- rotational  states  of 
molecules  (Shore,  1973;  1975;  Davis  and  Heller,  1979). 

For  very  high,  diffuse  (but  bound)  levels,  basis  sets  tend 
to  be  quite  large  and  there  is  parameter  adjustment  which 
must  be  done,  since  one  wishes  to  describe  both  very 
spatially  localized,  low-lying  levels  as  well  as  these 
highest  states  in  the  same  basis;  usually  an  improvement 
in  the  description  of  the  highest  levels  comes  at  the 
expense  of  the  lower  eigenvalues.  For  scattering  states  a 
still  larger  basis  set  would  be  required.  There  would 
also  be  a "termination"  error  associated  with  the  fact 
that  the  scattering  states  do  not  obey  boundary  conditions 
appropriate  to  bound  states;  we  must  put  our  particles  in 
a "box"  to  do  the  basis  set  calculation.  In  this 
connection  one  could  employ  a scatter ing- like  basis  (e.g., 
oscillatory  functions)  and  now  the  termination  error  would 
apply  to  the  bound  levels  rather  than  the  positive  energy 
states.  The  problems  discussed  here  are  technical  rather 
than  theoretical,  however,  and  it  re  mains  to  be 
demonstrated  whether  one  can  implement  this  methodology  to 
obtain  "bound  states  in  the  continuum"  of  nuclear  motion. 
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Returning  to  our  discussion  of  the  structure  in  the 
electron  energy  spectrum,  the  explanation  of  peak  (e)  is 
readily  obtained  with  the  use  of  Figure  4.  Here  we  depict 
the  initial  and  final  state  effective  potentials  for  the 
dominant  partial  wave  (L  = 24)  and  note  that  the  final 
channel  energy  lies  below  the  barrier  maximum;  this  is  a 
predissociating  state,  and  one  could  in  principle  obtain 
its  width  and  lifetime.  States  of  this  kind  have  long 
been  observed  in  diatomic  molecule  absorption  spectra; 
unlike  sharp  bound-bound  lines,  the  decay  of  the 
predissociating  state  results  in  a broadened  spectral 
line,  reflecting  the  finite  lifetime  of  the  excited 
species.  It  is  not  surprising  that  there  is  a marked 
enhancement  in  the  PEED  under  these  conditions;  the  longer 
the  particles  remain  bound,  the  greater  the  likelihood  of 
electron  emission.  Before  leaving  the  discussion  of 
electron  spectra  it  is  difficult  tc  resist  the  following 
speculation:  rather  than  "fold"  our  computed  results  with 
(unknown)  distribution  functions  (to  account  for 
ncn-monoenergetic  atomic  beams)  and  "acceptance"  functions 
(to  account  for  detector  resolution),  it  is  instead  hoped 
that  as  the  experimental  work  in  this  area  becomes  more 
refined,  observation  of  quantum  effects  such  as  those 
discussed  in  this  section  will  be  able  to  provide  a more 
stringent  test  of  the  theory. 

At  this  point  we  would  like  to  elaborate  on  a 
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possible  reaction  "mechanism"  for  Penning  ionization.  A 
study  of  molecular  orbital  correlation  diagrams  (Bellum 
and  Micha,  1977)  pictured  the  ionization  as  a molecular 
Auger-type  process,  going  through  a quasimolecular 
intermediate  from  which  an  electron  is  ejected.  In  the 
specific  case  cf  He*(1s2s,3S)  + Ar,  a simple 

interpretation  has  a "hole"  in  He*  going  over  into  a hole 
in  the  guasimo  lecule;  an  electron  from  a higher  fT 

molecular  orbital  (with  predominantly  Ar  3p  character) 
"falls"  into  the  hole  and  the  outermost  orbital  (singly 
occupied  and  strongly  He*  2s  in  character)  is  ionized 
resulting  in  a molecular  ion. 

Alternatively,  we  might  ask  whether  electron  emission 
actually  occurs  from  the  quasi  molecule  AB*,  or  from 

another  state  more  closely  resembling  (AB*,e~)  . This 
question  arises  from  the  fact  that  we  have  rationalized 
the  PEED  structures  in  terms  cf  g uasimolcular  states  in 

the  final,  molecular  ion  channel,  not  in  the  initial 

channel;  electron  emission  is  favored  when  HeAr*  is 
temporarily  trapped  in  a centrifugal  barrier.  Initial 
channel  resonances  could  certainly  be  used  to  explain 
enhancements  to  the  cross  section  from  certain  partial 
waves  as  the  colliding  particles  would  be  held  by  a 
centrifugal  barrier  long  enough  for  electron  emission  (or 
exchange  plus  emission)  to  occur.  To  test  this  hypothesis 
we  have  examined  the  initial  channel  for  centrifugal 
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barriers  and  have  found  that  1)  for  L > 18  the  effective 
potentials  are  purely  repulsive,  and  2)  for  L < 13  the 
largest  barrier  height  was  twenty  times  less  than  the 
collision  energy  (65  meV)  and  that  these  maxima  occurred 
at  internuclear  separations  greater  than  about  11  bohr. 
Thus  there  is  little  likelihood  of  initial  channel 
resonances  of  this  type  affecting  the  cross  section  as  in 
Figure  3.  Initial  channel  resonance  behavior  would  be 
more  likely  to  contribute  to  PEED  structure  and  larger  PI 
cross  sections  at  much  lower  collision  energies  (in  light 
of  remark  2)  above)  and  it  would  be  interesting  to  extend 
the  work  presented  here  to  much  lower  energies  to  test 
this  hypothesis.  A rise  in  the  ionization  cross  section 
at  lower  energies  would  be  consistent  with  the 
semie  mpirica 1 work  of  Olson  (1972)  and  which  has  already 
been  seen  in  PI  cf  Kr  by  Ne*  (Neynaber  and  Magnuson,  1976) 
and  of  Ar  by  Ne*  (*icha  et  al.,  1971;  Tang  et  al.  , 1972). 

Consequently  we  are  left  with  the  notion  of  final 
channel  resonances  as  an  explanation  for  our  PEED  results 
and  we  conclude  from  this  that  a HeAr*  -like  state  must  be 
formed  .prior  to  electron  ejection;  the  system  may  evolve 
as 

He*  ♦ Ar  -->  HeAr*  (guasi-molecule  formation) 

HeAr*  -->  (HeAr*  , e~)  (Auger  process) 

(HeAr*  , e~)  -->  HeAr+  * e~  (electron  emission) 

HeAr*  ♦ e~  — > He  + Ar+  * e~  (break-up) 
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with  final  channel  resonance  contributions  in  the  third 
step.  The  final  channel  maximum  occurs  away  from  the 
potential  minimum  and  we  may  think  of  HeAr*  in  this  state 
as  already  approaching  its  asymptotic  limit  of  He  ♦ Ar*. 
Although  Rydberg-like  states  were  not  included  in  our 
electronic  expansion  basis  which  gave  rise  to  the  coupled 
channels  equations  it  seems  possible  that  the  transition 
may  proceed  through  a state  in  which  an  electron  is 
loosely  bound  to  a HeAr*  ground  state  core.  He  then  see 
a dramatic  effect  of  nuclear  motion  upon  ionization 
probability:  our  initial  state  undergoes  an  Auger-type 
transition  to  an  intermediate  resembling  an  electron 
loosely  bound  to  HeAr+,  and  through  continued 
interaction  facilitated  by  rotational  trapping — the 

likelihood  of  ionization  is  increased. 

At  this  point  it  is  desirable  to  exami ne--by  way  of 
numerical  studies  the  validity  of  the  decoupling 
approximation  of  Chapter  I;  e.g.  , how  do  the  the  results 
of  the  35  (decoupled)  2-state  calculations  compare  with 
those  of  the  exact  matrix  problem  of  dimension  36  ? The 
motivation  for  this  study  was  made  clear  in  the  first 
chapter;  recall  that  with  the  neglect  of  the 
continuum-continuum  couplings,  the  potential  matrix  0 
blocked  in  such  a way  as  to  make  decoupling  a physically 
appealing  approximation.  There  are  several  considerations 
which  mandate  caution  in  proceeding  to  solve  the 
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higher-order  matrix  differential  equation.  One  factor  is 
the  linear  dependence  problem  already  discussed;  there  is 
a strong  likelihood  that  non-classical  regions  will  impede 
efficient  numerical  solution.  We  could  suppress  certain 
solutions  as  has  been  already  suggested,  but  in  the  event 
of  propagation  through  centrifugal  barriers  this  is  not 
justified.  A more  practical  problem  is  the  anticipated 
expense  of  such  a study.  Besults  of  Allison  (1970) 
suggest  that  for  the  deVogelaere  algorithm  the  time  goes 
as  approximately  N2*s  where  N is  the  number  of  coupled 
equations  (here,  the  number  of  discretization  basis 
functions  plus  one).  With  this  figure  our  anticipated  CFU 
time,  per  collision  energy,  is  roughly  40  hours;  this 
study,  much  less  a PI  energy-dependence  study,  is  clearly 
impractical  using  this  technique.  Lester  (1971)  indicates 
that  the  deVogelaere  algorithm  is  to  be  preferred  to 
Numerov  and  RK  methods  based  upon  consideration  of 
stability,  storage  requirements  and  machine  time. 
Although  the  piecewise-analytic  approach  of  Gordon  (1971) 
was  not  included  in  the  analysis,  that  method  is  claimed 
to  have  an  N3  dependence;  unlike  numerical  methods,  the 
step  size  is  not  related  directly  to  the  collision  energy 
(or  deBroglie  wavelength)  so  that  fewer  steps  may  be 
required.  It  seems  unlikely,  however,  that  there  will  be 
adequate  savings  to  make  even  this  method  desirable  for 
our  applications.  In  the  following  section  we  discuss  the 
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i npleieentaticn  of  a rather  recently-developed  method  for 
the  solution  of  a set  of  coupled  first  order  ODE's  which 
will  have  practical  consequences  for  our  study  of  Penning 
ionization. 


Solution  of  a Matrix  Riccati  Equation 
In  the  preceding  section  we  demonstrated  the  need  for 
a new  computational  approach  if  we  wished  to  assess  the 
validity  of  the  decoupling  approximation  of  Chapter  I.  As 
it  is  only  asymptotic  information  (e.g.,  K matrix)  which 
is  of  importance  to  us,  any  equivalent  formulation  of  the 
problem  capable  of  giving  us  the  desired  information  is 
acceptable.  One  such  approach  is  given  here.  Transform 
the  matrix  SchtOdinger  equation  for  the  wavefunction 
(suppressing  partial  wave  indices) 

D X = 0 X (3.1) 

by  use  of 


y = X'  X-i  (3.2) 

(where  prime  denoted  differentiation  with  respect  to  R) 
into  a matrix  equation  for  the  logarithmic  derivative 

y * ♦ W + y*  = 0 (3.3) 


where 

W = (2J1/X2)  (El  - V - Z}  . (3.4) 

Here  V is  the  effective  potential  matrix  (including 
centrifugal  potenial) , Z a diagonal  matrix  of  channel 
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threshhold  energies  and  1 the  unit  matrix.  The  boundary 
condition  on  X at  the  origin  gives  for  y ( 0)  a diagonal 
matrix  with  infinite  elements;  owing  to  the  symmetry  of  V, 
y (R  >0)  is  symmetric.  Usual  numerical  techniques  for 
solving  (3.3)  cannot  be  applied  because  y diverges  for 
certain  values  of  the  argument  (whenever  X has  a node).  A 
derivation  of  an  algorithm  for  the  solution  of  (3.3)  which 
circumvents  this  problem  and  which  is  based  upon  a 
limiting  form  of  the  ’’amplitude-density"  equations  of 
Johnson  and  Secrest  (1966;  1968)  may  be  found  in  Appendix 
B.  The  first  published  application  of  the  method  to 
scattering  problems  (Johnson,  1973)  demonstrated  its 
advantages  over  methods  for  solving  the  second  order 
matrix  differential  equation  in  regard  to  selection  of 
starting  points,  propagation  through  classically  forbidden 
regions  and  computation  time.  A more  recent  extension  of 
the  method  has  been  to  the  bound-state  solutions  of  the 
one-dimensional  Schrbdinger  equation  (Johnson,  1977),  an 
extreme  test  of  stability  in  propagation  through 
non-clasical  regions. 

Returning  to  the  asymptotic  form  of  X (1.26)  we  can 
write  for  K : 


K = - (yti-N’}-*  (yJ-J'}  (3.5) 

where  N and  J are  given  by  (1.27,1.28)  and  (1.35,1.36)  for 
open  and  closed  channels,  respectively.  Symmetry  of  y 
results  in  symmetry  of  K,  and  K may  be  again  partitioned 
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into  submatrices  (1.37).  In  the  following  section  we 
report  on  implementation  of  this  method  for  a study  of  the 
Penning  ionization  of  Ar  by  He*(1s2s,3S)  at  thermal 
energies. 

III.  3 Application  to  PI  of  Ar  by  He*(1s2s.3Sl 
at  Thermal  Energies 

A program  LOGDERIV  was  written  to  solve  the  Riccati 
transform  (3.3)  of  the  matrix  SchrOdinger  equation  for 
Penninq  ionization  in  atomic  systems.  Before  discussing 
applications  to  physical  problems,  it  is  appropriate  to 
relate  certain  observations  on  the  numerical  method.  As 
already  mentioned,  y(0)  is  a diagonal  matrix  with  infinite 
matrix  elements.  Johnson  (1973)  suggests  replacing  this 
by 

Y(0)  = cl  (3.6) 

where  c is  a large  number  (e.g.,  10*°) . We  examined  the 
asymptotic  solution  y for  a number  of  choices  of  c ranging 
from  1040  to  10-*°  and  in  all  cases  the  solutions  were 
identical  to  15  decimal  digits;  hence,  the  choice  of  y(0) 
is  quite  arbitrary,  and  this  is  a reflection  of  the  fact 
that  for  the  second-order  ODI  there  was  considerable 
flexibility  in  specif ication  of  the  first  derivatives. 
The  location  of  the  origin  was  varied  inside  the  classical 
turning  point  and  the  finding  was  the  same,  and  we 
conclude  (as  did  Johnson  (1973))  that  the  algorithm  is 
highly  insensitive  to  starting  conditions.  Another  point 
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we  note  is  that.,  from  the  results  of  Appendix  C,  the 
algorithm  requires  a matrix  inversion  at  every  integration 
step.  This  suggests  that  the  matrix  inversion  will  be  the 
time  consuming  operation  in  the  propagation.  For  the 
matrix  inverter  we  used  a routine  due  to  Beebe  (1978) 
which  is  based  upon  the  successive  orthonormalization 
procedure  of  LOwdin  (1968)  and  which  we  verified  through 
numerical  studies  to  be  an  N3  method.  Timings  were  made 
on  LCGDEBIV  using  this  rcutine  and,  based  upon  studies  of 
up  to  50  coupled  equations,  a very  clear  N3  dependence  was 
evidenced  here  as  well. 

Returning  to  the  present  application,  the  code  was 
tested  against  various  scattering  problems  in  the 
literature,  and  in  all  cases  the  results  were 
satisfactory.  For  elastic  scattering  we  again  compared 
with  the  phase  shift  results  of  Bernstein  (1960) ; for 
multichannel  problems  with  only  open  channels  we  compared 
with  H (or  K)  matrix  elements  from  an  exactly  soluble 
two-state  model  (Harris  and  Michels,  1969)  and  with 
transition  probabilities  for  He+  *■  Ne  inelastic 
scattering  (Johnson,  1973)  using  potentials  and  couplings 
of  Olson  and  Smith  (1971)  . The  inclusion  of  closed 
channel  boundary  conditions  was  compared  with  partial  wave 
cross  sections  of  Matese  and  Henry  (1972)  for  a model 
problem  which  has  some  of  the  features  of  low-energy 
electron-hydrogen  scattering. 
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For  PI  we  first  studied  the  convergence  of  the 
transition  probability  as  a function  of  discretization 
basis  size  for  the  "exact"  matrix  problem  (1.18).  The 
basis  was  that  of  Section  1.3  and  we  spanned  a fixed 
interval  of  the  electronic  continuum  with  an  increasing 
number  of  basis  functions.  in  Figure  5 we  present  a 
typical  convergence  study  for  PI  of  Ar  by  He*(1s2sr3S)  at 
a collision  energy  of  65  meV;  the  results  for  L = 0 
demonstrate  that  convergence  is  obtained  with  twenty  or  so 
basis  functions,  and  similar  behavior  was  noted  for  other 
partial  waves  and  at  other  energies.  Thus,  for  all 
further  results  to  be  reported  in  this  section,  the 
benchmark  was  a matrix  problem  of  dimension  20,  and  this 
will  be  compared  with  the  decoupling  results  of  19 
two-state  calculations. 

We  wish  to  examine  the  validity  of  decoupling  as  a 
function  of  two  dynamical  variables  in  the  theory,  namely 
the  heavy  particle  angular  momentum  partial  wave  L and  the 
collision  energy  E.  In  Figure  6 the  L = 0 partial  wave 
contribution  to  the  PEED  is  shown  for  a collision  energy 
of  65  meV  and  it  is  seen  that  the  agreement  between  exact 
and  decoupled  results  is  quite  good.  At  the  same  energy 
results  are  shown  for  L = 24  in  Figure  7,  and  in  Figure  8 
are  L = 0 results  at  125  meV.  Our  general  conclusion  from 
these  and  other  results  is  that  the  decoupling 
approximation  is  highly  satisfactory  and  that  its  accuracy 
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exhibits  no  dependence  on  collision  energy  (at  least  over 
the  range  of  interest)  or  on  heavy  particle  angular 
momentum  partial  wave.  The  results  in  Figures  6-3  are  not 
the  whole  story,  however,  as  one  must,  at  each  collision 
energy,  sum  over  L to  obtain  a PEED  (1.31)  and  integrate 
over  6 to  give  the  Penning  ionization  cross  section.  in 
Figure  9 we  show  a PEED  at  65  meV  and  discover,  in  keeping 
with  the  partial  wave  findings,  good  agreement  between 
exact  and  decoupled  results.  Although  the  continuum  mesh 
was  larger  here  than  in  the  study  in  Section  III.1,  the 
structure  in  Figure  9 again  arises  from  a dominant  partial 
wave  (L  = 24).  That  such  structure  was  evident  when  using 
an  entirely  different  computational  approach  lends  support 
to  the  interpretation  given  such  peaks  in  Section  III.1; 
they  are  decidedly  not  artifacts  of  the  calculation.  It 
is  highly  instructive  to  point  out  that  the  computations 
of  Figure  9 reguired  about  124  minutes  CPU  for  the 
20-state  case  and  about  8.7  minutes  for  the  decoupled 
calculations  on  an  Amdahl  470/V6  computer;  this  is  roughly 
8 seconds  versus  0.5  seconds  per  partial  wave  and  per  6. 
This  clearly  indicates  that  "benchmark"  calculations  are 
obtainable  at  a not  excessive  cost  and  that,  a cross 
section  energy-dependence  study  is  a feasible  project  to 
underta  ke. 

In  this  regard,  our  last  effort  in  the  present 
ccupled  channels  work  on  PI  is  to  study  the  Penning 


78 


ionization  cross  section  over  a range  of  collision 
energies  in  the  thermal  energy  range  (here  from  26  - 95 

meV)  where  a variety  of  experimental  work  is  available 
(Illenberger  and  Niehaus,  1975;  Pesnelle  et  al.,  1975; 
Wang  et  al.  , 1976;  Woodard  et  al.  , 1978);  our  results  are 
presented  in  Figure  10.  We  see  that  at  higher  energies  the 
rise  in  the  experimental  result  begins  to  decrease.  This 
may  be  understood  cn  rather  general  grounds  by  considering 
such  factors  as  the  duration  of  collision,  the  distance  of 
closest  approach  of  the  particles,  and  the  strength  of  a 
coupling  potential.  At  low  energies  the  particles  may 
spend  appreciable  time  in  ac  '•interaction  region"  but  at 
such  a distance  that  the  coupling  (or  transition)  strength 
is  quite  small.  As  energy  increases,  the  particles 
interact  for  a shorter  period  of  time  but  the  increased 
coupling  strength  at  shorter  distance  drives  the  reaction. 
At  still  higher  energies  (not  displayed  on  this  figure) 
the  particles  approach  closer  together  and  the  coupling  is 
stronger,  but  the  decreased  collision  time  is  the  dominant 
factor  and  the  cross  section  decreases. 

While  the  agreement  between  theory  and  experiment  is 
good  at  low  energies  it  clearly  deteriorates  at  the  higher 
energy  end.  There  are  several  factors  which  could 
influence  this  behavior:  1)  the  coupling  potential 

(Bellum  and  Micha,  1978)  was  parametrized  based  upon 
consideration  of  autoionization  widths  at  65  meV,  so  that 
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as  one  gets  too  far  away  fro®  this  energy  (higher  or 
lower)  it  is  not  likely  to  be  a realistic  coupling.  2)  We 
have  neglected  to  include  associative  ionization  channels 
in  our  coupled  channels  theory;  the  effect  of  this  process 
can  be  seen  in  the  experimental  results  of  Pesnelle  et  al. 
(1975) , the  AI  cross  section  is  small  at  low  energies  (20 
ineV)  , rises  to  a maximum  at  about  100  meV  and  falls  off 
again.  To  the  extent  that  PI  and  AI  are  competing 
processes  it  may  be  reasoned  that  the  inclusion  of  AI  in 
the  theory  might  tend  to  bring  the  entire  curve  down 
somewhat,  with  the  shift  being  greatest  at  the  high  energy 
end  of  Figure  10.  3)  The  reliability  of  the  potentials 

themselves  may  be  questioned.  This  aspect  of  the  problem 
has  been  discussed  earlier  and  we  reiterate  here  that  our 
emphasis  has  been  on  a detailed  study  of  the  dynamics  of 
the  heavy  particles  in  collisional  ionization,  with 
more-or-less  reasonable  interatomic  potentials  and 

couplings.  It  is  not  intended  that  the  potentials  employed 
here  should  be  capable  of  reproducing  a wide  variety  of 
experimental  results;  in  fact,  in  sem iclassical  and 

optical-model  studies  it  is  evident  that  high  and  low 
energy  data  cannot  be  reproduced  by  one  set  of  potential 
parameters  (e.g.  Woodard  (1978)).  Where  there  has  been  an 
attempt  to  "invert"  differential  elastic  scattering  cross 
section  data  to  obtain  the  relevant  potentials,  there  is 
considerable  discrepancy  between  the 


resulting  curves 
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(Brutschy  et  al.,  1976;  Wang  et  al.,  1976;  Altpeter  et 
al. , 1977;  Siska,  1979b).  it  is  hoped  that  the  discussion 
in  Section  II. 1 of  electronic  states  appropriate  to 
chemi-ionization  processes  will  serve  to  guide  theoretical 
efforts  in  this  regard. 
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I MEN:  McIntosh  Plot  System  - Mexico  16/03/80  23.10.59 
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Table  1 

Partial  Wave  Contributions  to  PEED  of  Figure  3 

L 


peak 

€ (a. u .) 

E (a.u.) 
f 

L 

% Contrib.  v (a.u. 

max 

a 

0.  1504388 

9. 5 (-4) 

29 

25 

9.219  (-4) 

b 

0.  1505388 

8-5  (-4) 

28 

27 

8.882  (-4) 

c 

0. 1506388 

7-5  (-4 ) 

27 

27 

7.495 (-4) 

d 

0. 1 507888 

6.  0 (-4) 

25 

24 

5 . 96  7 (-4) 

e 

0.  1508888 

5. 0 (-4) 

24 

56 

5 . 27  3 (-4) 

f 

0.1509883 

4.  0 (-4) 

20 

13 

2.94  3 (-4) 

19 

18 

2.472  (-4) 

18 

16 

2.046 (-4) 

g 

0. 1512388 

1 . 5 (-4) 

16 

27 

1.336  (-4) 
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CHAPTER  IV 

SOLUTION  OF  AN  INHOMOGENEOUS  DIFFERENTIAL  EQUATION 
FOR  PENNING  AND  ASSOCIATIVE  IONIZATION 

IV.  1 Approximations  to  a Complex,  Non-Local  Kernel 
We  have  derived  an  equation  for  initial  channel  heavy 
particle  motion  governed  by  a potential  which  is  complex, 
non-local  and  energy  dependent.  This  wavefunction  is 
needed  to  compute  various  scattering  cross  sections  of 
Section  III. 2.  Before  specifying  a method  for  solution  of 
the  inhomogeneous  second-order  ordinary  differential 
equation  (2.58),  we  will  in  this  section  examine  the 
different  contributions  to  the  kernel  to  see  whether 
simplifications  can  be  made.  In  the  next  sub-section  we 
review  some  traditional  approximations  to  the  complex 
kernel;  the  validity  of  these  approximations  will  be 
investigated  in  the  following  sub-section. 

IV. 2.1  Assumptions  Leading  to  a Local,  Imaginary 
Potential  Kernel 

Let  us  first  focus  our  attention  on  the  "electronic" 
contribution  to  the  kernel  from  PI  and  AI;  we  had  in 
(2.57)  the  terms 

p.v.  Jd€  V (R)  { E— € - E 1 } ~ 1 V (F')  - 

d€  €d 

-inV  ( R)  V (R • ) + 

dE'  ' E ' *d 
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p.v.  jde  V (R)  {E-G-E  }~1V  (B*)  - 

dG  L Gd 

-inV  (B)  V { R • ) 

dE"'  E"'d  (4.1) 

n 

with  E*  * = E-E*  and  E,,,;=E-E  for  PI  and  AI,  respectively. 

L 

In  the  first  term  E'  is  a continuous  variable 
corresponding  to  the  allowed  final  heavy  particle  states 

n 

of  PI;  in  the  third  term  E corresponds  to  a bound 

L 

vibrational  and  rotational  level  for  AI.  Two 

approximations  can  be  introduced  at  this  point  to  reduce 

the  complex  and  non-local  potential  QHPG+PHQ  to  an 

P 

imaginary  one  which  is  local.  The  first  step  is  to  assume 

(Bieniek,  1978;  Nakamura,  1972)  that  the  coupling  matrix 

element  V (R)  is  essentially  independent  of  G (the 
dG 

validity  of  this  assumption  will  be  examined  more  closely 

PI  AI 

in  the  next  sub-section);  the  terms  and  4 of  (2.62) 

L L 

are  zero  and  we  are  left  with  an  imaginary,  non-local  but 
separable  kernel 

K { E ; B , R * ) = -in  y (E  ; B,R  ' ) = -in* 

L L 


{ 5 d F»  U (E',P)V  (R)  V (R’)H  (£',R') 

I dE* ' E* ' d L 
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a 

+ 2(n)  v (E  ,R)  V * 

L L dE"»R 


n * 

V (R')f  (E  , R ' ) } 

E ' ' ' d L L (4.2) 

As  a further  approximation,  we  note  that  the  kernel 

contains  products  of  the  form  fV;  the  coupling  V decays 

exponentially  with  increasing  R and  the  functions  V 

oscillate  in  the  potential  well  region  and  then  either 

decay  (AI)  or  oscillate  (PI)  asymptotically.  This 

suggests  that  r|V  is  peaked  in  6 for  a value  of  e around 

the  Fr anck-Ccndon  value  € , with  the  result  that  V can  be 

v 


factored  out  of  the  sum  and  integral  in  (4.2)  . Using  the 

completeness  relation  for  the  set  of  heavy  particle 

functions  allows  us  to  write 

loc  r 

K ( E;  P , R ' ) = -inV  (R)  V (R  ' ) o (R- R ' ) 

L d€v  6vd  (4.3) 

where  we  choose  for  6 , using  electronic  adia baticity , 

v 


€ = e (R)  = V (R)  - V (R) 

V V d + (4.4) 

which  is  just  the  vertical  energy  difference  between  the 

resonance  and  ion  electronic  energies.  With  this  result 

equation  (2.58)  now  reads 

D X (E,R)  = -(2M/JTZ)|V  (R)  | 2X  (E,  R) 

I*  L d€v  L (4.5) 


or 


ii  o.i 
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{d*/dR2  - (2 M/Jf*)  (V  (E)  - iy  (R)  /2) 

d 

+ 2M  E/Jf2 } X (ErR)  = 0 
L 


(4.6) 


with 


y (R)  = 2tt  | V (R)  I 2 

d€w 


(4.7) 

electronically 


Although  originally  starting  from  an 
nonadiabatic  description,  the  result  (4.6)  has  reduced  (to 
within  neglect  of  the  real  part  of  the  complex  level 
shift)  to  the  result  of  Nakamura  (1969),  Fujii  et  al. 
(1970)  and  Miller  (1970)  where  an  adiabatic  description 
was  assumed  from  the  start. 


IV. 1.2  Real  Part  and  Non-Local  Nature  of  the 
otential  Kernel 

In  light  of  the  approximations  just  discussed,  it  is 
instructive  to  examine  the  contribution  of  the  p. v. 
integrals  to  the  kernel  and  thus  to  the  real  part  of  the 
optical  potential.  Quadrature  rules  for  p. v.  integrals 
have  been  developed  by  Elliott  and  Paget  (1979)  assuming 
the  functional  dependence  of  the  integrand  upon  the 
integration  variable  is  known.  However,  using  ab-initio 
electronic  wavef unctions  to  compute  coupling  matrix 
elements  and  autoionization  widths  is  a fairly  difficult 
task,  and  the  6-dependence  of  the  coupling  matrix  element 
is  generally  not  explicit.  One  approach  would  be  to  fit 
the  product  of  couplings  to  a function  of  6 for  which  the 
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p.v.  integrals 

can 

be 

evaluated  analytically. 

In 

the 

present  work 

we 

can 

obtain  an  estimate  of 

their 

contribution  in 

a 

manner  which  suffers  from  none 

of 

the 

above  drawbacks. 

We 

have 

already  seen  (Appendix  A) 

how 

the 

coupling  matrix  element  can  be  given  a simple  6-dependence 
based  upon  a Cculcmb  representation  of  the  continuum 


orbital;  for  this  case  the  p. v.  integrals  may  be  done 
analytically,  and  the  real  part  of  the  PI  contribution  to 
(4.1)  becomes 

PI 

A (R,R',E,E‘)  = (2/H)  (2 M/X2)"/2 
L 

rad  rad 

V (R)  V (R  ' ) { (E")-^ 

d€  6d  (4.8) 

ln|[  (€•+  (E")  )/(€•-  (E")  31-in  (E"  )~/2} 

A I n 

where  E''=E-E'  for  PI;  a second  term  A (R,F',E,E  ) is  of 

L L 

n 

the  same  form  but  where  E''=E-E  ; 6'  is  the  maximum 

L 

possible  value  of  6,  determined  by  energy  conservation, 
rad 

V (R)  is  the  radial  (R-dependent)  part  of  the  coupling, 
d€ 

given  in  Appendix  A.  Evaluating  (4.8)  and  its  AI 

counterpart  over  a range  of  E'  appropriate  to  PI  and  for 
n 

values  of  E corresponding  to  the  L=0  vibrational  states 
L 

of  the  ion,  the  real  part  was  found  to  be  of  the  same 
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magnitude  as  the  imaginary  part.  Furthermore,  over  this 
range  of  final  channel  energies  the  real  part  decreased 
rocnotonica lly  but  the  overall  change  was  less  than  1%. 
Whether  the  G-dependence  of  our  coupling  is  an  accurate 
one  would  need  to  be  demonstrated  by  actually  computing 
couplings  using  electronic  wa vef unctions  such  as  those 
suggested  in  Chapter  II.  In  this  connection,  Hickman  et 
al.  (1977)  observed  that  for  HeH*,  a 10%  increase 
(decrease)  in  k ( = €2/2)  causes  roughly  a 20%  decrease 
(increase)  in  the  autoionization  width.  For  our 

parametrizaton  of  the  coupling,  the  same  change  in  k 
results  in  about  a 10%  change  in  a computed  width. 
Although  a strict  comparison  is  of  course  unrealistic,  our 
results  indicate  an  S-dependence  of  roughly  the  same 
magnitude,  and  we  tentatively  conclude  that  evaluation  of 
p.v.  integrals  with  realistic  coupling  matrix  elements 
will  yield  a weak  final  channel  energy  dependence  similar 
to  that  observed  here. 


n 

Since  dependence  of  the  A on  E'  and  E was  observed 

L L 


to  be  weak,  the  A can  be  taken  outside  the  integral  over 

L 


E'  and  the  sum  over 
relation  for  the 
and  combining  terms 
p.v.  terms  to  the  po 


n in  (4.2).  Osing  the 
heavy  particle  functions 
we  obtain,  for  the  contri 
tential,  the  result 


completeness 
for  PI  and  AI 
bution  of  the 
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pv 

V (R,E, E ) = A ( R , R , E , E ) 

L v L v (4.9) 

where,  consistent  with  the  locality  approximation,  a 

plausible  choice  for  E is  that  obtained  by  assuming  a 

v 

Franc k-Ccndon  like  transition  for  the  nuclei  (i.e.,  one  in 
which  the  local  relative  kinetic  energy  of  the  nuclei  is 
conserved)  : 

E - V (R)  = E (P)  - V (R)  . 

d v + (4.10) 

Before  proceeding  to  discuss  the  imaginary  part  of  the 

kernel,  we  will  from  this  point  on  assume  that,  in  the 

general  case,  the  p.v.  integrals  are  obtainable  and  that  a 

local  approximation  to  the  result  has  been  added  to  the 

initial  channel  potential  V (R) . The  relative  magnitude  of 

d 

pv 

V as  compared  with  V has  no  bearing  upon  future 
L d 

developments.  The  key  step  in  the  above  treatment  was  the 

recognition  that  a local  approximation  to  the  p.v. 

results  seems  to  be  quite  reasonable.  If  we  were  not  able 

to  make  this  local  approximation,  we  would  have  to  carry 

along  a non-local,  real  contribution  to  the  kernel. 

Implications  of  this  are  discussed  in  the  next  section. 

Turning  our  attention  to  y'(E;R,Rl),  the  imaginary 

L 

part  of  the  kernel  composed  of  PI  and  AI  contributions, 
the  local  approximation  may  also  be  studied.  This  is  a 
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much  more  severe  approximation;  whereas  the  real  part 

(local  or  non- local)  of  the  kernel  is  essentially  a 

perturbation  to  V , the  imaginary  part  is  responsible  for 

d 


decay  by  electron  emission.  A local  approximation  here 
disallows  the  possibility  of  any  effect  of  the  relative 
motion  upon  the  ionization.  In  light  of  our  interpretation 
of  the  PEED  structure  in  Chapter  III,  there  is  support  for 
the  idea  that,  at  least  for  PI,  electron  ejection  may  have 
strong  dependence  upon  the  relative  nuclear  motion.  We 
may  ask  whether  the  same  can  be  said  for  AI  where  the 
heavy  particles  remain  bound  after  ionization  and  the 
electron  energy  spectrum  is  discrete.  To  help  understand 


AI 

the  AI  process  we  would  like  to  calculate  y (E;R,R»)  given 

L 

from  (2.57)  by 


AI  n 

Y (E  ; R,  R ' ) = -in  2 (n)  f (E  ,R) 
L L L 


n * 

*V  (R)  ip  (E  , R • ) V (R  • ) 

dE'"  L L E’"d  (4.11) 

n 

with  E,,,=E-E  and  where  the  sum  runs  over  all  vibrational 
L 

levels  of  the  ion,  in  our  case  HeAr*.  To  this  end  a 
program  VIBROT  was  written  to  determine  bound  state 
solutions  to  the  radial  SchrtJdinger  equation.  Following 
LeRoy  and  Bernstein  (1968),  Cashion  (1  963)  and  Cooley 
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(1961),  at  a trial  energy  a solution  was  integrated  inward 
from  some  "asymptotic"  point  until  the  outermost  maximum 
is  located;  another  solution  is  propagated  outward  from  an 
interior  point  to  the  point  cf  outermost  maximum.  Both 
solutions  are  scaled  to  unity  at  the  matching  point.  A 
correction  to  the  energy,  based  upon  the  difference  in 
slopes  at  the  match  point,  is  calculated  and  the  process 
is  repeated  until  convergence  is  achieved;  the  solution  is 
then  normalized.  The  vibrational  quantum  number  is 
identified  by  counting  the  nodes  in  the  solution.  In  our 
application  the  Numerov  method  was  chosen  as  the 

integrator.  The  program  was  successfully  tested  on  the 
bound  states  of  a Morse  oscillator,  whose  eigenvalues  are 
known  analytically  for  L=0. 

The  49  bound  vibration-rotation  levels  of  He Ar+  are 

given  in  Table  2.  In  this  connection  we  note  that  7 

quasibound  levels  were  found;  these  are  summarized  in 

Table  3.  Quasibound  states  here  are  those  which  are 

temporarily  trapped  in  a centrifugal  barrier  but  whose 

energy  lies  above  the  threshhold  fcr  dissociation.  These 

states  have  a lifetime  t and  a width  (energy  uncertainty) 

r 

y related  by 
r 

t = K/y  . 
r r 


For  these  resonance  levels  the 


outermost 


(4.12) 

boundary 
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condition  was  that  the  wavefunction  behave  as  an  Airy 
function  of  the  second  kind  (Abramowitz  and  Stegun,  1972) 
at  the  third  turning  point.  This  Airy  function 
exponentially  increases  in  extreme  non-classical  regions, 
is  easy  to  calculate  for  small  arguments  (such  as  near 
turning  points),  and  was  found  to  result  in  more  accurate 
eigenvalues  than  other  choices  for  outer  boundary 
conditions  (LeRoy  and  Bernstein,  1971).  The  semiclassical 
approximation  to  t based  upon  the  vibrational  period  and 


r 

barrier  tunneling  probability  given  in  LeRoy  and  Bernstein 


(1971) 

was  used 

to  calculate  the 

results  in  Table 

3. 

These 

results 

are  presented  here 

because  in  addition 

to 

"pure  " 

AI  there 

may  be  long-lived 

molecular-ion  states 

whose 

detectio  n 

depends  upon  their  lifetime  and 

on 

experimental  detection  time.  At  a typical  detection  time 
of  10-6  seconds  (Hickman  and  Morgner,  1976)  only  the  n=  1 , 
L=14  level  is  likely  to  be  "observed"  as  dissociated 
products;  most  of  the  quasibound  states  can  thus  be 
treated  as  true  AI.  The  existence  of  such  states  suggests 
that  the  sum  over  bound  states  appearing  in  this  theory  be 
modified  to  include  quasibound,  positive-energy  states  as 
well.  For  these  energies  the  electron  energy  spectrum 
corresponding  to  these  levels  will  overlap  with  the  PEED; 
there  will  be  discrete  lines  embedded  in  the  electron 
energy  distribution  at  the  corresponding  energies. 


In  Figures  11  and  12  we  plot  the  vibrational 
wavef unctions  for  HeAr*  for  L=0,  n=0-5.  What  is 
surprising  is  the  spatial  extent  of  the  levels  n=3,4,5; 
recall  that  the  position  of  the  potential  minimum  is  about 
6 bohr.  It  is  of  concern  to  us  that  the  long-range 
wavef unctions  will  contribute  to  the  imaginary  part  of  the 
kernel  and  thus  to  the  effective  potential.  However,  we 
see  that  the  kernel  is  a product  of  bound  wavef unctions 
with  potential  coupling  matrix  elements;  the  latter  have 
an  exponentially  decreasing  R-dependence  and  we  expect 
that  any  contributions  from  the  long-range  wavef unc tions 
will  be  minimized. 

Feturning  to  our  discussion  of  the  imaginary  part  of 
the  kernel  for  AI,  we  have  plotted  in  Figures  13-18,  at  a 
collision  energy  of  65  meV,  the  individual  terms  (to 
within  a factor  -in)  in  the  sum  over  n in  (4.2);  since  the 
kernel  is  symmetric  in  R and  R’  the  "hidden"  part  of  the 
figures  can  be  readily  visualized.  For  n=0  the  kernel  is 
everywhere  positive  since  the  vibrational  wavefunction  is 
nodeless.  For  n=1-5  only  two  undulations  are  observed,  and 
these  are  at  distances  in  the  region  of  the  potential 
well;  the  exponentially  decaying  coupling  damps  out  any 
oscillatory  behavior  at  large  R.  In  Figure  19  the  "sum 
over  states"  kernel  is  shown;  it  is  dominated  largely  by 
the  n=0,1  vibrational  levels  since  only  their 


wavef unctions  have  maxima  in  this 


region. 


Based  upon 


Ill 


these  results  it  is  not  likely  that  a local  approximation 
to  the  individual  terms  in  the  AI  kernel  is  warranted. 

In  conclusion,  our  experience  with  Penning  ionization 
(Chapter  III)  suggests  that  nuclear  motion  does  play  a 
role  in  the  electron  emission  process;  this  emission  is 
caused  by  the  imaginary  part  of  the  kernel,  and  any  local 
approximation  to  this  quantity  will  undoubtedly  be  a poor- 
one.  For  AI  an  analysis  of  the  terms  in  the  imaginary 
part  of  the  kernel  suggest  that  a local  approximation  will 
be  likewise  inappropriate. 

We  will  new  direct  our  attention  to  a method  of 

solution  of  an  inhomogeneous  differential  equation  for 

Penning  and  associative  ionization  with  a complex, 

separable  kernel.  Recalling  that  the  p.v.  contribution  to 

the  kernel  may  be  approximated  in  a local  manner  and 

combined  with  the  initial  channel  potential  V , 

d 

simplifications  arising  from  an  imaginary,  non-local 
kernel  will  also  be  discussed. 

IV. 2 Method  of  Solution 

Before  proceeding  with  a method  of  solution  of  (2.58) 
it  is  of  interest  to  point  out  a simplification  which 
results  from  our  formulation  of  the  effective  Hamiltonian 
problem  in  Chapter  II.  We  had  ottained  an  inhomogeneous 
second-order  ordinary  differential  equation  (ODE)  by 
making  a spectral  decomposition  of  the  operator  g*(t=)  in 


1 12 


terms  of  heavy  particle  functions  for  PI  and  AI;  this 
choice  ensures  a separable  potential.  Had  we  used  the 

closed  form  for  the  radial  Green's  function  (Micha,  1976) 
we  would  have  likewise  obtained  a separable  potential  but 
for  an  inhomogeneous  integro-dif ferential  equation  (IDE) 
for  the  heavy  particle  wavef unction.  A non-iterative 
solution  of  the  latter  requires  that  one  solve  a 

homogeneous  IDE  plus  one  (or  more)  particular 
inhomogeneous  IDE.  In  the  present  approach  we  will  see 
that  we  will  be  required  tc  solve  a homogeneous 
second-order  ODE  plus  a collection  of  inhomogeneous  ODE's; 
the  construction  of  the  kernel  appears  simpler  with  the 
present  formulation. 

Our  method  of  solution  of  the  non-local  potential 
problem  for  complex  potentials  is  in  the  spirit  of  the 
work  of  Gordon  (1970)  and  Percival  (as  communicated  by 
Marriott  (1958)).  We  first  review  the  method  for  real 
kernels  to  establish  notation. 

Consider  the  radial  equation 

D ( R)  f (R)  = J (R  , R ) dR*  K(R,R')f(R') 

1 2 (4.13) 

where 

D(R)  = (d2/dR2  ♦ (2M/tf2)  (E-V  (R)  ) - L(L  + 1)/B2}  (4.14) 

and  where  we  have  suppressed  energy  and  partial  wave 
arguments  for  clarity.  For  a separable  kernel 

K ( R ,R  ' ) = £(i=1,N)  V (R)  W ( R 1 ) 

i i 


(4.15) 
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then 


Df  = £(i)  V <W  |f> 

i i (4.16) 

where  < | > denotes  an  integral  over  R'.  The  quantities 

<W  |f>  are  constant,  independent  of  R and  so  we  have  an 
i 

inhomogeneous  equation  with  a fixed  functional  form  on  the 
right  hand  side.  First  we  solve 


and  the  set 


Oh  = 0 


(4.17) 


DI 


i 


V 

i 


(4. 18) 


for  i=1  to  N.  Now  the  most  general  solution  f can  be 

written  as  a particular  solution,  which  is  a linear 

combination  of  the  solutions  I , plus  any  homogeneous 

i 


solution  (such  as  h)  : 

f = h + Z.(i)  C I . 

i i (4.  19) 

It  is  readily  shewn  by  substitution  of  (4.19)  into  (4.16) 

that  the  coefficients  {C  } satisfy  the  set  of  simultaneous 

i 

linear  equations  (in  matrix  notation) 


(1  - M )C  = M 

1 2 (4.20) 

where  1 is  the  unit  matrix  of  dimension  N,  C is  the  vector 

of  coefficients,  the  matrix  H is  given  by 

1 


(M  ) = <W  |I  > 

1 ij  i 1 


(4.21) 
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and  M is  the  vector 
2 

(M  ) = <W  |h>  - 

2 i i (4.22) 

The  solution  h is  easily  obtained:  it  is  just  the 

one-dimensional  analog  of  the  solution  to  the  matrix 

SchrOdinger  equation  discussed  in  chapters  I and  III.  For 

the  inhomogeneous  solutions  the  reader  is  directed  to 

Gordon  (1970)  where  a piecewise-a na lytic  method  for 

obtaining  solutions  such  as  I is  discussed. 

i 

For  the  case  of  a complex,  separable  kernel  we  write 

again 

D(R)f(R)  = j (0rinf)  dfi*  K ( R , R 1 ) f (R  ' ) (4.23) 

but  now  for  K we  choose  the  expression 

K ( R , R ' ) = Z.(i=  1 , 2N)  V (P)  W (R») 

i i (4.24) 

where  V and  W are  real  for  i= 1 to  N;  for  i=N+1  to  2N  we 
i i 

choose  V real  and  make  W imaginary.  As  before  we  solve 
i i 

Dh  = 0 (4.25) 

and 


DI  = V 

i i (4.26) 

for  i=1  to  2N.  By  specifying  W rather  than  V to  be 

i i 

imaginary  for  i=N+1  to  2N  those  solutions  I can  all  be 

i 

real  if  we  choose  K-matrix  boundary  conditions  as  in 


Chapter  I.  Me  write  again 


f = h + Z.(i=1,2N)  C I 

i i 


(4.27) 


and,  as  before,  the  {C  } satisfy 


i 


(1  - H )C  = H 


2 


(4. 28) 


where  now  the  matrices  1 and  H (4.21)  and  the  vector  H 

1 

(4.22)  are  of  dimension  2N.  Since  some  elements  in  M can 


now  be  imaginary,  the  solutions  C are  in  general  complex. 
A substantial  simplification  has  been  achieved  in  this 
formulation:  by  making  W be  the  imaginary  part  of  the 


kernel  for  i=M  + 1 to  2N  we  have  deferred  operations 

involving  complex  arithmetic  until  the  last  step  (4.28). 

Had  we  chosen  to  make  V (i=N<- 1 to  2N)  imaginary,  the 

i 

inhomogeneous  solutions  to  (4.26)  would  be  complex. 

We  now  proceed  to  apply  the  above  method  to  the 
homogeneous  ODE  discussed  in  the  previous  section. 
Inspection  of  the  AI  part  of  the  kernel  (2.57)  allows  us 
to  write  (suppressing  energy  arguments  and  partial  wave 
indices  for  clarity) 


2 


l 


AI 

R = Z(n=1,N)  u v * Z(m  = 1,N)  w x 


(4.29) 


with  the  identifications 


(4.30a) 


1 16 


n 

u = V (E  , R)  * (R-dependent  part  of  p.v.) 
n 

n 

v = |(E  ,R')  * (R '-dependent  part  of  p.v.) 
n (4.30b) 

m 

w = r (E-E  ,R)  V (R) 

m dE'"  (4.30c) 

m 

x = f (E-E  ,R  ' ) (-in)  V (R») 

m E'"d  (4.30d) 


{ E * * • = E-E  ) or,  combining  in  a single  sum 


A I 


where 


for  p=1  to  N and 


for  p=N+ 1 to  2N. 


PI 


K = Z.(p=1,2N)  0 V 

P P 


0 = u 

P P 


V = v 
P P 


U = w 
p p-N 

V = x 
p p-N 


(4.31) 

(4.  32a) 
(4 . 32  fc) 

(4.  32c) 
(4. 32d) 


For  K we  make  a quadrature  approximation  to  the 
integral  over  E':  using  N'  points, 

PI  c 

K = JdE*  K(E')  = 2.(r=1,N‘)  K (E  )z 

r r r 
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= Z(r=1  ,N')  K'  (E  ) 

r r (4.33) 

where  (z  } are  quadrature  weights  which  in  the  last  line 


have  been  absorbed  into  K* 

r 


the  form 


PI 

Now  K has  been  brought  to 


PI 

K = 2(r=1,»')  a b ♦ Z(s=1,N')c  d 

r r s s 

where 

a = If  (E  , R) * (R-de penden t part  of  p.v.) 
r r 

b = i|(E  , R •)*  (R  ' -dependent  part  of  p.v.) 
r r 

C = $(E-E  r R ) V (R) 
s s dE* ' 


d = ? (E-E  ,F*  ) ( in)  V (R  • ) 
s s E'  1 d 

(E'  • = e-E  ) or 
s 

PI 

K = Z(q= 1 , 2N 1 ) A B 

q q 

with  the  identification 

A = a 

q q 

B = b 

q q 

for  q=1  to  N'  and 

A = c 
q q-  N ' 


(4.34) 

(4.35a) 
(4.35b) 
(4. 35c) 
(4. 35d) 


(4.36) 


(4. 37a) 
(4.  37b) 

(4.37c) 
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B = d 
q q~N ' 


(4.  37d) 


for  q = N • +1  to  2 N ' . 

We  now  have  the  result 


or 


K = Z(p=1,2N)  U V + Z(q=1,2N')  A B 

P p q q 


K = Z(t=1,T)  F G 
t t 


where 


and 


for  t=1  to  2N,  and 


T = 2N  + 2N ' 


F = U 

t t 

G = V 

t t 


F = A 

t t 

G = B 

t t 


(4.38) 

(4.39) 

(4.40) 

(4. 41a) 
(4.41b) 

(4.41c) 


(4. 4 Id) 

for  t=2N  * 1 to  T.  With  so  toe  effort  we  have  been  able  to 
put  the  kernel  for  AI  and  PI  into  a form  (4.39)  such  as 
(4.24)  which  we  can  solve  in  exactly  the  same  manner. 
Obtain  the  solutions  to 

D h = 0 (4.42) 

DI  = F 

t t (4.43) 

for  t=1  to  T.  Then  the  solution  X of  the  inhomogeneous 
differential  equation 
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DX  = 5 KX  (4. 44) 

= Z(t=1,T)  F <G  |X> 
t t 

is  given  as 

X = h + Z(t=1,T)  I Y 

t t (4.45) 

and  the  coefficients  (Y  } are  the  solution  to  the  matrix 

t 

eguat ion 

( 1 - H ) Y = H 

1 2 (4.46) 

where  the  unit  matrix,  M and  M are  of  dimension  T. 

1 2 

Specifically, 


(H  ) 

= <G  | I > 

1 

ij 

i j 

(4.47) 

(H  ) 

= 

<G  | h>  . 

2 

i 

i 

(4.48) 

The 

expression 

(4 

.4  5) 

for 

the  solution  of  (4.44) 

is  the 

main 

result  of 

this 

section. 

The  functions  G are 

either 

i 

real  or  imaginary,  and  the  solutions  I are  real;  thus  the 

i 

elements  of  M and  a are  either  real  or  imaginary,  but 
1 2 

not  complex.  As  we  found  above,  we  are  able  to  defer  any 
arithmetic  operations  involving  complex  numbers  until  the 
very  last  step  (4.46)  , the  solution  of  the  set  of  linear 
equations  which  determine  the  coefficients.  Although  this 
last  point  is  not  critical  to  the  mathematical 

development,  it  is  certainly  of  practical  importance  when 
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doing  machine  computation:  manipulations  directly  in  terms 
of  complex  quantities  exact  an  appreciable  expenditure  of 
time  due  to  the  increased  number  of  arithmetic  operations 
required. 

In  light  of  the  discussion  in  the  previous  section 

regarding  approximations  to  the  full  kernel,  in  those 

cases  where  we  may  either  neglect  the  principal-value 

contributions  or  where  a local  approximation  way  be  made 

to  the  real  part  of  the  kernel,  the  following 

simplifications  are  noted:  the  quantities  u , v , a and 

n n r 

b are  absent,  and  the  dimensionality  of  our  matrix 
r 

problem  is  reduced  by  a factor  of  two  (now  T = N + N ' ) - 

a and  M are  composed  of  imaginary  elements,  and  in  this 
1 2 

case  the  complex  coefficients  Y can  be  obtained  by 

i 

solving  a real  set  of  simultaneous  linear  equations. 

This  concludes  our  discussion  of  a method  for  the  (in 

principle)  exact  solution  of  an  inhomogeneous  differential 

equation  arising  from  a theory  of  ionizing  collisions 

between  atoms.  Once  we  have  constructed  the  homogeneous 

solution  h and  the  particular  solutions  I , the  M arrays 

i 

can  be  calculated  and  the  desired  scattering  solution 
obtained  algebraically.  In  the  last  section  we  will 
summarize  the  principal  results  of  the  present  work  and 
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discuss  problems  where  fuxther  studies  are  warranted. 


IV. 3 Discussion 

Our  study  of  ionizing  atomic  collisions  at  low 
energies  has  caused  us  to  focus  on  both  electronic  and 
heavy-particle  dynamic  aspects.  The  electronic  portion  of 
the  problem  concerns  the  description  of  the  appropriate 
electronic  wave  functions,  potential  energy  curves  and 
coupling  matrix  elements.  Discussion  was  given  of  the 
construction  of  a continuum  orbital  using  the  expansion 
method  of  Harris  (1967)  which  has  enjoyed  widespread 
application  for  electron-atom  scattering.  Continuum 

electronic  wavef unctions  are  then  built  up  from  ground 
state  molecular  ion  wave f unc tions  and  these  continuum 
orbitals.  We  also  proposed  hew  to  obtain  molecular 
resonance  or  "hole"  wavef  unctions;  a resonance 

wavefunction  which  is  explicitly  orthogonal  to  the 
continuum  electronic  functions  was  obtained  by  projection 
operator  technigues.  Alternatively,  an  approximate 
resonance  wavefunction  was  obtained  from  a ground  state 
function  by  the  appropriate  single- part icle  excitation; 
this  function  was  then  or thogonalized  to  the  continuum 
electronic  functions  to  give  an  orthogonal  basis  of 
electronic  wavef unct ions  suitable  for  our  scattering 


problem. 
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»e  made  a detailed  numerical  investigation  of  Penning 
ionization  of  Ar  by  He* ( 1s2s, 3S)  using  a coupled-channels 
formalism  (Bellum  and  Micha,  1978)  for  the  heavy  particle 
wavef unctions . Detailed  analysis  of  the  energy 
distribution  of  the  emitted  electron  suggested  a profound 
effect  of  nuclear  motion  on  the  ionization  probability. 
Appreciable  enhancement  of  this  energy-differential  cross 
section  was  found  to  occur  at  energies  corresponding  to  at 
least  one  heavy  particle  guasibound  state  and  perhaps 
several  "discrete  states  in  the  continuum"  of  nuclear 
motion.  This  led  us  to  propose  that  the  initial 
electronic  state  undergoes  a transition  to  an  intermediate 
state  resembling  an  electron  loosely  bound  to  an  ionic 
core;  it  is  from  this  intermediate  that  the  electron  may 
then  be  ejected.  A method  was  described  whereby  these 
"discrete"  continuum  heavy  particle  states  could  be 
obtained  by  applying  basis  set  expansion  techniques  to  the 
scattering  states  as  well  as  to  the  bound  states  of 
nuclear  motion. 

The  validity  of  approximating  the  full  matrix 
SchrOdinger  equation  in  a coupled  channels  theory  by  a 
collection  of  two-state  scattering  equations  was 
demonstrated  by  exploring  alternative  computational 
techniques;  here  we  chose  to  solve  the  Riccati  transform 
of  the  matrix  SchrOdinger  equation  using  a method  which 
afforded  a significant  reduction  in  computer  time.  The 
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decoupling  approximation  was  tcund  to  be  independent  of 
energy  and  heavy  particle  angular  momentum  partial  wave 
over  the  range  of  energy  of  interest.  Penning  ionization 
cross  sections  were  obtained  and  found  to  compare 
favorably  with  experimental  results  at  energies  from 
20-100  meV. 

In  Chapter  II  an  effective  Hamiltonian  theory  was 
developed  for  collisional  ionization  along  the  lines  of 
Feshbach's  unified  theory  of  nuclear  reactions  (1958; 
1962)  . Using  projection  operator  techniques  we  derived  an 
effective  potential  for  the  colliding  particles  which  was 
complex,  non-local  and  energy-dependent.  In  Chapter  IV  a 
detailed  analysis  was  made  of  the  Penning  and  associative 
ionization  contributions  to  the  kernel  of  a inhomogeneous, 
second-order  ordinary  differential  equation.  He 
demonstrated  that  although  the  real  part  of  the  kernel  may 
by  simplified  by  a local  approximation,  it  was 
inappropriate  to  similarly  treat  the  imaginary  part.  He 
also  discussed  formation  of  long-lived  nuclear  motion 
states  in  the  final  channel  and  suggested  that  their 
lifetimes  be  determined  so  that  they  may  be  correctly 
classified  as  PI  or  AI ; this  classification  of  course 
depends  upon  an  experimental  detection  time.  A method  was 
developed  for  the  solution  of  such  inhomogeneous  equations 
with  complex  kernels;  the  method  is  general,  and  we 
discussed  its  i nplementa tion  for  PI  and  AI  and  pointed  out 
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simplifications  resulting  from  the  case  of  a purely 
imaginary  kernel.  A particular  advantage  of  the 
formulation  defers  any  operations  involving  complex 
quantities  until  the  final  step,  which  would  have 
practical  consequences  for  an  actual  computation. 

There  are  still  many  interesting  aspects  of 
collisional  ionization  which  merit  further  attention.  For 
example,  extension  of  either  the  coupled  channels  cr 
effective  Hamiltonian  theory  to  include  higher  angular 
momentum  partial  waves  for  the  ejected  electron  is 
desirable.  Although  limiting  the  treatment  here  to  the 
1=0  electron  partial  wave  is  advantageous  from  the  point 
of  view  of  simplifying  the  study  of  nuclear  dynamics, 
higher  partial  waves  of  course  must  contribute.  A further 
improvement  of  either  theory  would  be  the  inclusion  of 
multiple  electronic  continua.  For  the  He*-Ar  system 
studied  here,  there  is  a 2n  molecular  ion  electronic  state 
which  may  also  be  formed.  This  feature  gives  rise  to 
interference  effects:  the  initial  discrete  state  is  now 
coupled  to  two  continuum  electronic  states,  and  there  are 
in  principle  couplings  between  the  two  continua  as  well  as 
within  each  continuum. 

From  the  computational  point  of  view,  we  mentioned  in 
Chapter  III  that,  in  addition  to  final  channel  heavy 
particle  resonance  effects,  initial  channel  quasibound 
states  may  modify  the  dynamics  at  lower  collision 
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energies.  It  would  be  of  interest  to  see  whether  the  cross 
section  rises  at  ouch  lower  energies  and  if  such  a rise 
could  be  explained  by  resonance  behavior.  A more  exciting 
prospect  is  the  implementation  of  the  method  for  solving 
inhomogeneous  equations  for  complex  kernels  as  detailed  in 


this  chapter. 

The  solution  requires  homogeneous 

a nd 

inhomogeneous 

solution  s 

as  input. 

and  these  may 

be 

obtained  by 

standard 

methods. 

Application  to 

an 

analytically  solvable  two-state  scattering  problem  would 
be  a first  step  toward  demonstrating  its  feasibility  and 
in  affording  comparisons  with  traditional  coupled-channels 
treatments. 
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Table  2 

Heftr(+)  Bound  Levels 

E ( 1/cm) 

L 

n=0  n=1  n=2 

n=3 

n=4 

0 

- 1.409  (2)  -5.048(  1)  -1.637(  1) 

-4. 182 

-5.  73  (-1) 

1 

- 1.400(2)  -4.978  (1)  -1.591  ( 1) 

-3. 913 

-4.608  (-1) 

2 

-1.  382  (2)  -4.838  ( 1)  - 1.  499(  1) 

-3.  386 

-2.533  (-1) 

3 

- 1.  356  (2)  -4.629  (1)  -1.363  (1) 

-2.670 

a 

4 

- 1.320  (2)  -4.  352  (1)  -1.  183  (1) 

-1.  652 

5 

- 1.276  (2)  -4.007  (1)  -9.  633 

-5.  400  ( 

-1) 

6 

-1.  223  (2)  -3.597  (1)  -7.  066 

7 

- 1.  162  (2)  - 3.  122  ( 1)  -4.  102 

8 

- 1.091  (2)  -2.586  (1)  - 1.  06  1 

9 

- 1.  012  (2)  - 1.992(1)  a 

10 

-9.250  (1)  - 1.343  (1) 

1 1 

-8.  293  (1)  -6.  456 

12 

-7.252  (1)  a 

13 

-6.128(1)  a 

14 

-4.924  (1)  a 

15 

-3.641  (1) 

16 

-2.286  (1) 

17 

-8. 448 

18 

a 

19 

a 

n=5 


-8.8  (-3) 


a = quasibound  level 
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Table  3 


HeAr{  + ) Quasibound  Levels 

n 

L 

E (1/cm) 

Width  (1/cm) 

Lifetime  | 

0 

18 

6.634 

3.5  (-20) 

1.5  (8) 

19 

2.240  (1) 

4.2  (-12) 

1.3 

1 

12 

9.  190  (-1 ) 

5.4  (-19) 

9.8  (6) 

13 

8.540 

2.  1 (-8) 

2.  5 (-4) 

14 

1.609  (1) 

2.6  (-6) 

2.0  (-6) 

2 

9 

2.  138 

4 . 0 (-8) 

1*3  (-4) 

4 

3 

2.30  (-3) 

2.4  (-12) 

2.2 
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APPENDIX  A 

INTERACTION  POTENTIALS  AND  COUPLINGS 
For  convenience  the  potentials  used  in  the  work  of 
Relyea  and  Micha  (1979)  and  Bellum  and  Micha  ( 1978)  are 
reproduced  below.  To  give  a simple  functional  form  to  the 
potentials  for  HeAr*  and  HeAr  + they  have  both  been 
represented  by  a Born-flayer  exponential  repulsion  plus  a 
long-range  van  der  Waals  attractive  tail.  They  can  be 
represented  by 

-R/a 

V (R)  = Ae  + f(R;RO,b,c)V  (R) 

LR  (A  . 1 ) 

where  f is  the  3-parameter  switching  function  of  Eckart 
(1930) 

(R-RO) /c 

f (R ; R° , b,c)  = (1  + e )-i* 

(R-E°) /c 

{1  + b (1  ♦ e )-»}  (A. 2) 

which  takes  on  the  values  (0,1)  at  the  endpoints  of  the 
interval  (0,°°  ) of  R,  and  which  exhibits  a local  maximum 
(for  b>0)  at 

R = RO  + c* ln[  (b+  1)/ (b-  1)  ] . 
max  (A. 3) 

The  switching  function  is  used  to  smoothly  join  the 

repulsive  and  attractive  parts  of  the  potential,  with 
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parameters  determined  by  reproducing  the  well-depth  and 

position  of  the  potential  minimum.  To  determine  R°  the 

value  of  R was  chosen  to  coincide  with  the  position  B 
max  m 

of  the  potential  minima.  The  Born-Mayer  parameter  a is 

deduced  from  the  MSX^  results  of  Bellum  and  Micha  (1977) 

on  HeAr*  and  HeAr*;  the  value  of  c is  set  equal  to  a.  To 

determine  A and  b we  use  the  conditions 

V (R  ) = -E 

m m (A. 4) 

V' (S  ) =0 

m 

where  S is  the  well  depth.  For  HeAr*  the  values  of  E 
m m 

and  R are  these  of  Chen  et  al.  (1974) ; those  relevant  to 
m 

HeAr4-  are  from  Nakamura  (1975).  Specifically,  we  obtain 
(energies  in  Hartrees,  B in  Eohr) 

-R/0. 9675 

V = 4.  39678e  + f (R  ; 9.4 48 8,  2.  8 1 505 , 0. 9675) 

HeAr* 

* (-2  26. 0/B6)  (A.  5) 

and 

-R/0. 532 

V = 44.  6527e  ♦ f (R  ; 5.  93 6,  6. 42 3 1 7 ,0. 532) 

H eA  r + 

* (-0.6904/R*-7.94/R*)  . (A. 6) 

The  coefficient  of  R-6  for  HeAr*  is  from  Bell  et  al. 
(1968);  for  HeAr*  we  use  the  R~*  coefficient  of  Dalgarno 
(1962)  and  the  R~6  coefficient  of  Nakamura  (1975). 
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The  coupling  matrix  element  V (S)  has  been 

d6 

parametrized  in  terms  of  an  autoionization  width 

y (R)  = 2n  I <♦  I H !<!>(€)  >|  2 

d el  (A. 7) 

which  is  often  given  the  exponential  representation 

-R/a° 

Y (R)  = A°e  (A.  8) 

in  experimental  and  semiempirical  analyses.  To  give  the 

model  coupling  an  explicit  6-dependence  (which  is  absent 

in  y of  (A. 8))  it  is  noted  that  the  normalization  constant 

for  a Coulomb  orbital  of  energy  6 is  (2/nk)  where 

k2=2m6/H2  (m=electron  mass) . We  write 

<t»  (€)  = (2/nk)  yH'  (A . 9 ) 

and  give  the  coupling  matrix  element  <4>  |H  | <J» • > an 

d el 

R-dependence  suggested  by  (A. 8).  Thus 

-P/a°  -R/a' 

A°e  = 2n  ( 2/nk)  (A  'e  )*  . (A.  10) 

Selection  of  A0  and  a0  values  from  semiempirical  analyses 

allows  us  to  obtain  values  for  A'  and  a'  and  the  coupling 

becomes 

-R/a' 

V (R)  = (2/nk ) I/2  A ' e 

dC  (A. 11) 

Tabulation  of  values  of  various  workers  may  be  found  in 

Bellum  and  Micha  (15  7 8)  . For  all  work  reported  here,  the 

values  A'  = 50.0  Hartrees  and  a'  = 0.7  Bohr  were  used. 
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APPENDIX  B 

TRANSITION  AN  PL IT  ODE  S T 

QQ 


T 

PC 


Recall  we  had  written 


(E  - QHQ  - QHEG  + PHQ) Of  + = 0 
P Q 


or 

(E  - H •)  Cf  ♦ = 0 
0 

where 

H*  = QHQ  * QHPG  +PHQ  = H + F* 
P QQ 

= I4>  > (T  ♦ V ) <♦  | + F+ 

d R d d 


= H°  ♦ V°  + V1 

with  the  choice  of  Q as  a projector  onto  the  closed 
electronic  channel, 

Q = 14  ><<D  I 
d d 

and  where  f is  the  total  wavefunction  (nuclear  and 

Q 

electronic)  when  the  initial  channel  is  Q.  We  also 

remember  that  H accounts  for  elastic  scattering  in  Q by 

QQ 

potential  V , and  F + accounts  for  transitions  into  the  P 
d 

subspace  (i.e.  , ionization)  via  interaction  PHQ, 
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propagation  in  P via  G ♦,  and  a final  interaction  QHP 

p 

which  accounts  for  recombination  in  P and  which  returns 
the  nuclei  to  their  initial  states  in  the  entrance  channel 
Q.  The  two- potential  formalism  (Rcdberg  and  Thaler,  1967) 
for  a potential  v = V°  + V1,  which  requires  eigenfunctions 
) a>  of  HQ,  | b>  of  H°  + V°,  and  | c>  of  H ' gives,  for 

transitions  from  channel  i to  channel  f, 

T = <a  | V 0 | b +>  ♦ 
fi  f i 

< b ~ | V i | c + > . 
f i 

In  our  case  we  make  the  identifications  (using  the 
sub-index  "q"  to  denote  functions  defined  only  within  the 
Q subspace)  : 


1 a>  = 

1?  °>  = 

1 >|u«> 

q 

d 

b+>  * 

||I+  > = 

|(j)  >|u  + > 

q d 

I c * > = | Q f ♦> 

Q 

= t<o  ><♦  ||)+  > 

d d Q 

= l«  > I X +> 
d C 

and  u+  and  u°  are  related  by  the  Lippmann-Schwinger 
equation 

u + = u o + G + V u + 

Q , 0 cl 


with  solution 
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u ♦ = u°  +G+V  u° 

Q d 

and 

g ♦ = {E  - t - v + ie)-i 
Q R d 

and 

G ♦ = (E  - T + i€)-»  . 

Q,0  a 

The  last  two  quantities  are  related  by 

G*=G  ♦ + G +VG+ 

Q Q# 0 Q , 0 d Q 

= G ++G+VG  + . 

Q, 0 Q d Q,0 

We  thus  obtain 

T = <V  o l <|>  >V  <4»  | V +> 

QQ  q d d d q 

+ <V  -|F+|Qf  ♦> 

q c 

= <u° | V Iu+>  + — | F* | QV  *>  * 

d q Q 

Recalling 

FQf  ♦ = QHPG  +PHQQI  ♦ = QHP|I  + 

Q P Q Q 

we  obtain,  alternatively, 

T = <u°|V  | u* > + <|  -IQHPIV  +>  . 

QQ  d q Q 

Note  that  the  first  term  on  the  right-hand  side  is  just 

the  transition  matrix  element  for  elastic  scattering  in  Q, 

uncoupled  to  P,  so  we  write 

T = T 0 + <f  ~ | F+ 1 Q|J  + > 

QQ  QQ  q Q 
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= T 0 ♦ <?  ~ | QH P | f ♦> 

QQ  g Q 

where 

T 0 = <u°  | V | u*>  . 

QQ  d 

The  distor ted-wave  Born  approximation  (DWB)  amounts 
to  replacement  of  eigenfunctions  of  H'  by  eigenfunctions 
of  H°  + VO  so 

Q V * — > Q I + = V * 

Q q q 

and  the  first  form  of  T becomes 

QQ 

DWB,  1 

T = T o + <|j  — |F  + ) |l  ♦> 

QQ  QQ  q q 

and  at  this  point  coupling  to  the  ionized  channel  is  still 

present  through  the  operator  F+. 

In  the  second  form  of  T using 

QQ 

Pf  ♦ = 9 * * QV  ♦ 

Q Q Q 

and  the  DWB  approximation  gives 

DWB, 2 

T = T 0 + <|  - J QH| V +> 

QQ  QQ  q Q 

- <?  -IQHQIP  + > • 

q q 

, the  ionization  amplitude,  we  make  use  of 
PQ 

Q(S  - H)Q|  = QHPf 
Q Q 

P (S  - H)  Pf  = PHQ |I 
Q Q 


For  T 
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and  the  corresponding  homogeneous  equations 

Q(E  - H)Q|  = 0 

q 

P(E  - H)  P|  = 0 
P 

to  write  the  channel  Hamiltonians 

H = QHQ  + QHP  = H°  + V 
Q 

and 

H = PHP  + PHQ  = HO ' * V*  . 

P 

Then  using 

T = <1  ”| V* | ¥ *> 

fi  f i 

gives 

T = <|  — | PHQ | 9 +>  - 
PQ  P Q 

The  DWB  approximation  here  gives  directly 

DWB 

T = <¥  -| PHQ | V ♦> 

PQ  p q 

so  again  we  need  not  determine  V or  its  P or  Q 

Q 

projections.  This  approximation  retains  states  fully 
distorted  by  elastic  scattering  in  Q and  P (f  ♦ and  | 

q p 

resp.)  and  is  superior  to  the  normal  Born  approximation  to 

T which  replaces  Qlf  f by  IJ  a free- motion  wavef  unction 
PQ  Q q 

(an  eigenfunction  of  H°)  which  is  not  distorted  by  elastic 

scattering  by  V in  Q space, 
d 


APPENDIX  C 

DERIVATION  OF  THE  LOG-DERIVATIVE  ALGORITHM 


Consider  the  radial  equation 

u"  (r)  + g (r)  *u  (r)  =0  (C.  1) 

where 

q (r)  = (2m/)*2)  ( E - V (r)  ) 

ef  f 

or 

q (r)  = k2  - o (r) 

and  where 

V (r)  = v (r)  + (K2/2  m)  1 (1  + 1 ) /r2 

eff 

For  p a positive,  real  number  write  (Johnson,  1973) 

u"  (r)  - p2  u ( r)  = -(  q (r ) + p2  )u(r)  (C.2) 

or 

u'»(r)  - p2u(r)  = - f(r)u(r)  (C . 3 ) 

where 

f(r)  = q(r)  ♦ p2  (C.  4) 

Consider  also  the  corresponding  homogeneous  equation 

w"  (r)  - p2w  (r)  = 0 

Then  the  solution  to  (C.1)  can  be  written  as 

pr  c -p(r-r') 

u(r)=e  + 1/(2p)  J (0,inf)  e f (r  1 ) u ( r • ) dr  • (C.5) 

Let 
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1 55 


F ( r)  = f (r)  *u  (r) 

then 

pc  f -p(r-r') 

u(r)=e  + 1/  (2p)  J^inf)  e F(r')dr' 

and 

lira  u (r)  = 

r-->i nf 

pr  -prc  pr* 

e + 1/(2p)e  J (0,in£)  e F(r')dr' 

Let 


K = 0 (0  ,i  nf ) e F (r)  dr 

Then 


pr  -pr 

lim  u (r)  = e + 1/(2p)e  K 

c — >inf 

Multiply  (C.5)  by  f (r)  ; then 

pr  c -p(r-r') 

F(c)=£(r)e  + 1/ (2p)j  (0,  inf)  e f(r)F(r')dr' 

Our  objective  is  to  obtain  (C.10),  with  K given  by 

and  F(r)  by  (C.11).  Perforin  the  integrals  in  (C 

(C - 11)  numerically,  using  N- 1 points  r ( j)  and 

c ( j)  . So 


pr  ( j) 

K = 2L(j=1,N-1)  e F (r  (j)  )c(j) 

1 1 


and 


pr  (i) 


(C.6) 

(C.7) 


(C.8) 

(C.9) 

(C. 10) 

(C.11) 
(C.9) 
.9)  and 
weights 

(C.  12) 


F (r(i))=f  (r{i))e 
1 


+ 1/(2p)  Z(j=1,N-1)  f (r (i) ) 
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-p  (r  {i ) -r  ( j)  ) 

*e  *F  (r  ( j)  ) c ( j) 

1 (C . 13) 

Now  repeat,  using  N points  and  weights: 

Pt  ( j) 

K = Z(j  = 1 ,N)  e F (r  (j)  ) c { j)  (C.  14) 

and 

pr  (i) 

F (r  (i)  ) =f  (r  (i)  ) e + 

-p  (r  (i)  -r  (j) ) 

1/(2p)Z(j  = 1,H)  f (r  (i)  )e  F(r(j))c(j)  (C . 15) 

Consider  (C.15)  at  r=r(N)  : 

F(r(N))  = f (r  (N) ) + 1 / ( 2p)  f ( r (N ) ) F (r  { N)  ) c { N)  * 

-p  (r  (N)  -r  (-j) ) 

1/(2P)Z(  j=1,N-1)  f (r  (N))  e F (r  ( j)  )c  (J)  (C.16) 

and  then 

F (r  (N) ) = (C.  17) 

pr  (N) 

( f (r  ( N)  ) e ♦ 1/  ( 2 p ) 2^j=1,N-1)  f(r(N))* 

-p  (r  (N)  -r  ( j)  1 

e F(r(-j))c(j)  ) * (1-1/  ( 2 p)  f(r(N)  )o(N))  1 


so 


pr  (N) 

F (r  ( N)  ) = -2p/c(N)e  A-1/c(N)A  Zj(j=1,N-1)  F (r  ( j) ) c ( j) 

-P(r  (Ni  -r  ( j) ) 

*e  (C.  18) 

where 

A=1/(2p)f  ( r(N))c(N)*(1/(2p)f(r(N))c(N)  -1)-*  (C.19) 

Thus  we  can  write 

pr(N)  -pr  (N) 

F (r  (N)  ) = -2p/c  (N)  e A - 1/c(N)Ae  b (C.20) 
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with  b to  be  determined. 


He  now  assert  that  for  i<N,  F(r(i))  is  related  to 

F (r(i))  by 
1 


F(r(i))  = F (r  (i) ) a 

1 (C  . 2 1 ) 

with  a to  be  determined.  Substitution  of  (C.20)  and 
(C . 2 1 ) into  {C.  15)  , 

pr  (i) 

F (r  (i)  ) a = f (r  (i)  ) e + 1/ (2  p)  t (r  { i)  ) Z(j=1,N-1)  ac(j) 

1 


-p(r  (i) -r  (j) ) -p(r  (i)-r(N) ) 

F {r  (j ) ) ♦ 1/(2p)  f (r  (i) ) e 
1 


pr(N)  - p r ( N ) 

{ -2p/c  ( N)  e A - 1/c  (N)  Ae  b)c(N)  (C.22) 

followed  by  substitution  of  (C.13)  into  the  left-hand  side 
of  (C.22)  leads  to 

pr  (i)  -p  (r  (i) -r  ( j)  ) 

(f(r(i))e  +1/(2p)Z.(j  = 1,N-1)  f(r(i))e  (C.23) 

pr  (i) 

*F  (r(j))c(i))a  = f (r  { i)  ) e + 1/ (2p)  f (r  (i)  ) 

1 


-p  (r  (i)  -r  (i) ) 

*2L(j  = 1,N-1)  e F (r  ( j)  ) ac  ( j)  + 1/  (2p)  *f  (r  (i)  ) 

1 


-p  (r  (i) -r  (N)  ) pr(N)  -pr(N) 

*e  (-2p/c  (N)  e A- 1 /c  ( N)  Ae  b)c(N) 

Simplifying,  we  get 

pr(i)  pr  (i)  -p  (r  (i) -r  (N) ) 

f (r  (i)  ) e a = f (r  (i)  ) e -f(r(i))e  * 

pr(N)  -p  (r  (i) -r  (N)  ) -pr(N) 

e A - 1/ (2 p)  f (r  (i)  ) e Abe  (C.24) 
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This  holds  for  all  r(i),  so  that  at  r(i)=r(N): 

-2pr  (N) 

a = 1 - A - 1/  {2p>  Abe  (C.  25) 

Now  we  obtain  b;  substitute  (C.21)  aud  (C.20)  into  (C.15) 
for  r (i)  =r  (N)  : 


pr(N)  -pr(N)  pr(N) 

-2p/c(M)e  A - 1/c  { N)  Ae  b=f(r(N))e  (C.26) 


-p  (r  (N)  -r  (j) ) 

+ V(2p)  £(j=1,N-1)  f (r  (N)  ) e F (r(j))ac(j) 

1 


pr(N)  -pr(N) 

♦ 1 / (2  p)  f { r (N ) ) (-2p/c{N)e  A-1/c(N)Ae  b)c(N) 

and  for  r(N)  > r ( j)  we  obtain 

pr(N)  -pr(N)  pc  (N) 

-2p/c(N)e  A - 1/c(N)Ae  b=f  (r  (N) ) e ♦ 


-pr(N)  pr(N) 

1/ (2p)  f (r  (N) ) e K a-f  (r  (N) ) (e  A ♦ 

1 


-pr  (N) 

1/  (2p)  e Ab)  (C.  27) 

Simplifying  (C.27): 


pr  (N)  -pr  (N) 

A (e  + 1/  (2p)  e b)  (-2p/c  (N) +f  (r  ( N)  ) = 

pr  (N)  -pr  (N) 

f (r  (N)  ) (e  >1/{2p)  e K a) 

1 (C.  28) 

Note  that 


-2p/c  (N)  + f (r  ( N)  ) = f (r(N))/A 

so  that 

b = K a 

1 

Now  express  K (C.  14)  in  terms  of  K (C.12): 

1 


(C . 29 ) 
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pr  (j) 

K=2T(j=1,S)  e F(r(j))c(j) 

pc(j)  pr(N) 

= Z(j=1,N-1)  e F (r  ( j)  ) ac  ( j) +e  F(r(N))c(N) 

1 


pr  W pr(N)  -pr(N) 

=K  a<-e  c (N)  (-2  p/c  ( N)  e A-1/c{H)e  Ab) 

1 


2pr  (N) 

=K  a - 2pe  A-Ab 


1 


or 


2pr  (N) 

K=b-2pe  A-Ab  (C 

Keeping  in  mine!  the  results  (C.25),  (C.29)  and  (C 

define 


-2pr  (N) 

b'=-1/(2p)e  b (C 

-2pr  (N) 

K =-  1 / (2  p)  e K 

2 (C 

-2pr  ( N—  1 ) 

K = -1  / (2p) e K 

3 1 


or 


2pr  (N-1) 

K =-2pe  K 

1 3 

2pr  (N)  -2ph 
~ - 2pe  e K 

3 (C 

where  h = r(M)  - r(N-1).  Then  (C . 25)  becomes 

a=1-A*Ab»  = 1 +A  (b ' - 1 ) (C 


. 30) 

.30)  , 

.31) 

. 32) 


33) 

34) 


(C.29)  becomes 
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-2  ph 

b ' =e  K a 

3 (C.  35) 

and  (C. 30)  becomes 

K =b  ' + A-Ab  ' = A + b ' (1-A) 

2 (C.36) 

Also  define 

- 2ph 

K =e  K 

4 3 (C.  37) 

so  that 

b • -K  a 

4 (C.  38) 

Substitute  (C.34)  into  (C.38)  to  obtain 

b'  = ( 1-K  A)  ~ 1 K (1-A) 

4 4 (C.  39) 

Substitution  cf  (C.39)  into  (C.36)  gives 

K =A«-  (1-K  A)-iR  (1-A)  2 
2 4 4 (C.  40 ) 

From  (C.10)  we  can  now  write,  at  r = r(N), 

pr  (N)  -pr(N)  2pr  ( N) 

u(r(N))=e  <-1/(2p)e  (“2pe  K ) 

2 


pr  ( N) 

=e  (1-K  ) 

2 


Differentiating 


(C.  10)  with  respect  to  r gives 


pr  -pr 

u ' (r)  =pe  +1/(2p)(-p)e  K 


pr  -pr 

= pe  +p(-1/(2p)e  K) 


(C.41) 


which  becomes,  at  r=r  (N ) 
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pr  (N) 

u ' (r  (N) ) =pe  ( 1 + K ) 

2 (C.42) 

Define  L{r)  , the  reciprocal  of  the  logarithmic 
derivative,  by 

L (r)  =u { r)  (u'  (r)  ) - * 


From  (C.  19)  , 


and  by  (C.36)  then 


Let 


= (1-K  ) (p  { 1 +K  ))-» 

2 2 


lim  A = 1 

p — > 0 


lim  K = 1 

p— >0  2 


Z=(1-K  ) / ( 2 p ) 

2 


(C.  43) 


(C.44) 


so 


Also  define  Z*  by 


lim  Z = L 

p — > 0 


K = 1-2pZ ' 
3 


t hen 


-Z=-{1-K  )/(2p) 

2 

= ( 1- A>  - (1  -K  A ) ~ 1 K ( 1 -A)  ? / (— 2p) 

2 2 


and 


(C.  45) 


(C.  4 6) 


(1-A)/  (-2p)=-1/(2p)  (1- (f(r(N))c(N))  (f  (r(M))c(N)-2p)-» 

= (f  (r  { N)  )c{N)-2p)~l 
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Now  expand  K in  powers  of  p,  keeping  terms  through  first 
4 


order : 


-2ph 

K = e K 
4 3 


= (1-2ph)  ( 1 -2pZ ' ) 
= 1-2p(h  + Z») 


Then 


1-K  A=  (f  (r  (N)c  (N)  -2p-  ( 1-2p  (h+Z' ) ) E (r  (N)  c (N)  ) 
4 


* (f  (r  (N)c  (H)-2p) -*  (C.4  7) 

and  after  straightforward  but  tedious  algebra  we  arrive  at 
Z=(h+Z'  + 2p(h+Z')/(f  (r(N)c(N))  (-1+2p/f  (r  (8)c  (N) 

* { 1-f  (r  (N)  c (N)  (h+Z')  ) (C.  4 8) 

and 

lim  Z = L = (h  + Z')  (1-q  (r  (N)  ) c (N)  (h+Z'))-1 

p — >0  (C.  49) 

and  new 


Z=  (h  + Z ' ) (1-g  (r  (N)  ) c (N)  (h  + Z'))-1 

Finally,  let 


y (i)  =z  (i)  -1 

and  we  obtain  the  result 


(C . 50) 


(C.  5 1 ) 


y (i)  = (hy ' + 1)  -1  y ' -q  (r  (i)  ) c (i)  (C.  52) 

where  y*  = y(i-1)  and  y (k)  =y  (r  (k)  ) . We  could  at  this  point 
take  c(i)  to  be  Simpson's  rule  weights,  i.e., 
c (i)  = h/3  at  y=y  (0)  and  y=y  (8)  , 

= 4h/3  at  y = y (1 ) , y (3)  , . . . ,y  (N-1 ) , 
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= 2h/3  at  y=y(2)  ,y  (4)  . ,y  ( N)  , 

where  N=even  integer.  Equation  (C.52)  is  the  log 
derivative  algorithm  presented  by  Johnson  (1973),  and 
additional  discussion  on  the  choice  of  c (i)  can  be  found 
in  Johnson  and  Secrest  (1966;  1968). 
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